Google 



This is a digital copy of a book that was preserved for generations on library shelves before it was carefully scanned by Google as part of a project 

to make the world's books discoverable online. 

It has survived long enough for the copyright to expire and the book to enter the public domain. A public domain book is one that was never subject 

to copyright or whose legal copyright term has expired. Whether a book is in the public domain may vary country to country. Public domain books 

are our gateways to the past, representing a wealth of history, culture and knowledge that's often difficult to discover. 

Marks, notations and other maiginalia present in the original volume will appear in this file - a reminder of this book's long journey from the 

publisher to a library and finally to you. 

Usage guidelines 

Google is proud to partner with libraries to digitize public domain materials and make them widely accessible. Public domain books belong to the 
public and we are merely their custodians. Nevertheless, this work is expensive, so in order to keep providing tliis resource, we liave taken steps to 
prevent abuse by commercial parties, including placing technical restrictions on automated querying. 
We also ask that you: 

+ Make non-commercial use of the files We designed Google Book Search for use by individuals, and we request that you use these files for 
personal, non-commercial purposes. 

+ Refrain fivm automated querying Do not send automated queries of any sort to Google's system: If you are conducting research on machine 
translation, optical character recognition or other areas where access to a large amount of text is helpful, please contact us. We encourage the 
use of public domain materials for these purposes and may be able to help. 

+ Maintain attributionTht GoogXt "watermark" you see on each file is essential for in forming people about this project and helping them find 
additional materials through Google Book Search. Please do not remove it. 

+ Keep it legal Whatever your use, remember that you are responsible for ensuring that what you are doing is legal. Do not assume that just 
because we believe a book is in the public domain for users in the United States, that the work is also in the public domain for users in other 
countries. Whether a book is still in copyright varies from country to country, and we can't offer guidance on whether any specific use of 
any specific book is allowed. Please do not assume that a book's appearance in Google Book Search means it can be used in any manner 
anywhere in the world. Copyright infringement liabili^ can be quite severe. 

About Google Book Search 

Google's mission is to organize the world's information and to make it universally accessible and useful. Google Book Search helps readers 
discover the world's books while helping authors and publishers reach new audiences. You can search through the full text of this book on the web 

at |http: //books .google .com/I 



A TREATISE 



ON 



THE STABILITY OF MOTION. 



■v* 



A TREATISE 



ON THE 



STABILITY OF A GIVEN STATE 

OF MOTION, 



PARTICULARLY STEADY MOTION. 



BEING THE ESSAY TO WHICH THE ADAMS PRIZE WAS ADJUDGED 
IN 1877, IN THE UNIVERSITY OF CAMBRIDGE. 



BY 



E. J. EOUTH, M.A., F.R.S., &c. 

LATE FELLOW OF ST PETEB'S COLLEOB, CAMBRIDGE, 
AND LATE EXAMINEE IN THE UNIVEBSITT OF LONDON. 



Uontfon : 
MACMILLAN AND CO. 

1877 

[All n'njhts reserved.] 



I?J 



s. . i,^ • 



■ «• . 









% 



^ . - . 



i 



> ■« 



I.*' 






> y ' 



■K- 



V 



* 






. r 

■ : 



(STatttbtiUge: 



PRINTEI^ BY 0. J. CLAT, U.A., 
AT THE UNITER8ITT PRESS. 



PREFACE. 



In March, 1875, the usual biennial notice was issued, giving the 
subjects for the Adams Prize to be adjudged in 1877. The 
following is the chief portion of the notice : 

The University having accepted a Fund raised by several members 
of St John^s College for the purpose of founding a Prize to be called the 
Adams Prize^ for the best essay on some subject of Pure Mathematics, 
Astronomy or other branch of Natural Philosophy, the Prize to be 
given once in two years, and to be open to the competition of all 
persons who have at any time been admitted to a degree in this 
XJuiversity — 

The Examiners give notice that the following is the subject of the 
Prize to be adjudged in 1877 : The Criterion of DTfnamical Stability, 

To illustrate the meaning of the question imagine a particle to slide 
down inside a smooth inclined cylinder along the lowest generating line, 
or to slide down outside along the highest generating line. In the 
former case a slight derangement of the motion would merely cause 
the particle to oscillate about the generating line, while in the latter 
case the particle would depart from the generating line altogether. 
The motion in the former case would be, in the sense of the question, 
stable, in the latter unstable. 

The criterion of the stability of the equilibrium of a system is, 
that its potential energy should be a minimum; what is desired is, a 

R. A. h 
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corresponding condition enabling us to decide when a dynamically pos- 
sible motion of a system is such, that if slightly deranged the motion 
shall continue to be only slightly departed from. 

The essays must be sent in to the Vice-Chancellor on or before the 
16th December 1876, &c., &c. 

S. G. PHEAR, rice-Chancellor. 

J. CHALLIS. 

G. G. STOKES. 

J. CLERK MAXWELL. 



The pressure of other engagements for some time prevented 
me from giving my attention to the subject. This essay was 
therefore almost entirely composed during the year 1876. It is 
now printed as it was sent in to the Examiners, the changes being 
merely verbal. Some few additions have been made where ex- 
planation appeared to be necessary, but all these have been 
marked by square brackets, so that they can be at once dis- 
tinguished from the original parts of the essay. 

In order to shorten the essay as much as possible many 
merely algebraic processes have been omitted and the results only 
are stated. It is hoped that this will add clearness as well as 
brevity to the reasoning, as the attention of the reader will not 
be called from the argument to follow a manipulation of symbols 
which may not present any novelty. 

Tte line of argument taken may be indicated in a general way 
as follows. Chapter i. begins with some definitions of the terms 
stable and steady motions. It is then pointed out that whether 
the forces which act on the system admit of a force-function or not, 
the stability of the motion, if steady, is indicated by the nature 
of the roots of a certain determinantal equation. The boundary 
between stability and instability being generally indicated by the 
presence of equal roots, a criterion is investigated to determine 
beforehand whether equal roots do or do not imply instability. 
This case being disposed of, the consideration of the determinantal 
equation is resumed. Two general methods are given by which, 
without solving the equation, it may be ascertained whether the 
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character of the roots imply stability or instability. These obcupy 
Chapters ii. and III. In the first method a derived equation is 
blade use of, and it is shown that a simple inspection of the signs 
of the coefficients of the several powers in these two equations will 
decide the question of stability. In the second method a certain 
easy process is found which if performed on the determinantal 
equation will lead to the criteria of stability. At the end of the 
third Chapter a geometrical interpretation is given to the argu- 
ment. 

In the fourth Chapter the forces which act on the system 
are supposed to have a force-function. The determinantal equa- 
tion is then much simplified. Several points are considered in 
this Chapter which are necessary to the argument, such as the 
proper method of choosing the steady co-ordinates (if there be 
any), the distinction between harmonic oscillation about steady 
motion and that about equilibrium, and the changes which must 
be made in the determinantal equation when the equations of 
Lagrange become inapplicable. A method of modifying the 
Lagrangian function is also given by which, in certain cases, 
the fundamental determinant may be reduced to one of fewer 
rows and columns. 

In the fifth Chapter a series of subsidiary determinants is 
formed, and it is shown that at least as many of the conditions 
V of stability are satisfied as there are variations of signs lost in the 
series in passing from one given state to another. It is also shown 
that this is equivalent to a maximum condition of the Lagrangian 
function. 

In the sixth Chapter the energy test of stability is considered. 
It is also shown that, when the motion is steady, this reduces 
to the same criterion as that indicated in Chap. V. 

In the seventh Chapter the question congftdered is whether 
the stability of a state of motion can really be determined by an 
examination of the terms of the. first order only. In some cases 
these are certainly sufficient, and an attempt is made to discrimi- 
nate between these cases and those in which the terms of the 
higher order ultimately alter the character of the motion. 

If the Hamiltonian characteristic and principal functions be 
given, the conditions of stability as regards sp^e only, or both 
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space and time may be deduced. But if these be not known as 
expressed in the Hamiltonian form, we may yet sometimes dis- 
tinguish between stability and instability if we can determine 
whether a certain integral ceases to be a minimum at some instant 
of the motion. This is the subject of the eighth chapter. 

As part of the third edition of my treatise on the Dynamics of 
Rigid Bodies was written at the same time as this essay, there are 
necessarily points of contact between the two works. Thus the 
subjects of the first part of the seventh chapter and of a portion 
of the sixth will be found discussed in the treatise on Dynamics. 
But as the objects of the two books are not the same, it will be 
found that in all these cases there are considerable differences in 
the modes of demonstration. 

EDWARD J. ROUTH. 



Peterhouse, 

August 14, 1877. 
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CHAPTEE I. 



Definitions of the terms small quantity, stable motion, steady motion. 
Arte. 1, 2. 

A system of bodies in steady motion is stable if the roots of a certain 
determinantal equation a/re stich that their real pa/rts a/re oM negative. 
Art 3. 

Effect of equal roots, and a test to determine whetlier equal roots do or do 
J^ not introduce terms which contain the time as afadLor. Arts. 4 — 7. 

Object of Chapters II, a/nd III. Art, 8. 

1. Let us suppose a dynamical system to be set in motion 
under any forces and to move in some known manner. If any 
small disturbance be given to the system, it may deviate only 
slightly from its known motion, or it may diverge further and 
further from it. Let 0, <f), &c. be the independent variables or 
co-ordinates which determine the position of the system, and let 
the known motion be given by ^ = ^^, <^ = <^o, &c. where 0^,^ (jy^, &c. 
are known functions of the time t To discover the disturbance of 
the system we put = 0q + x, ^ = ^o + y> &^- These quantities 
X, y, &c. are in the first instance very small because the disturb- 
ance is small. The quantities x, y, z, &c. are said to be small 
when it is possible to choose some quantity numerically greater 
than all of them, which is such that its square can be neglected. 
This quantity may be called the standard of reference for small 
quantities. 

If, after the disturbance, the co-ordinates x, y, z, &c. remain 
always small, the undisturbed motion is said to be stable; if, on the 
other hand, anyone of the co-ordinates become large, the motion is 
called unstable. 

It is clear that the same motion may be stable for one kind of 
disturbance and unstable for another. But it is usual to suppose 
the disturbance general, so that if the motion can be made un- 
stable by any kind of disturbance (provided it be small) it is said 
to be unstable. On the other hand, it will be called stable only 
when it is stable for all kinds of small disturbances. 

R. A. 1 



2 DEFINITIONS. [CHAP. 

2. To determine whether x, y, z^ &c. remain small, we must 
substitute for ^, ^, &c. in the equations of motion their values 
^0 + ^> <^o + y> ^c* Assuming that x, y, &c. remain small, we may 
neglect their squares, and thus the resulting equations will be 

linear in a?, y, Zy &c. The coefficients of a?, ^ , -^ » y» 7^ » 7/^' ^^• 

in these equations may be either constants or functions of the 
time. In the former case the undisturbed motion is said to be 
steady for these co-ordinates, in the latter unsteady. In the case 
of a steady motion a?, y, «, &c. are all functions of the time 
which has elapsed since the disturbance and of certain constants 
of integration which are determined by the initial values of 

doD dxi 

^y ~T±>y>~jiy ^^' ^® ^^y therefore define a steady motion to be 

such that the same change of motion follows from the same initial 
disturbance at whatever instant the disturbance is communicated 
to the system. 

If all the coefficients in the equations to find aj, y, z are con- 
stant, they may be made to contain ^ by a change of co-ordinates. 
Thus we may write for a?, y, «, &c. 

a? = af + /817 + . . . 

= &C. 

where a, /8, &c. are any functions of t we please. Conversely, when 
the coefficients are functions of ty we may sometimes make the 
coefficients constant by a proper change of co-ordinates. But this 
cannot always be done. If there are n co-ordinates, we have r? 
arbitrary functions a, /8, &c. at our disposal. In each of the n 
linear equations of motion we may have three terms for each 
co-ordinate, and thus we have (3 w — 1) w coefficients to make con- 
stants. We have therefore in general too many equations to 
satisfy. The proper method of choosing the co-ordinates of refer- 
ence will be considered in a future chapter. 

3. Let us suppose a dynamical system to be making small 
oscillations under the action of any forces which may, or may not, 
possess a force function and to be subject to any resistances which 
vary as the velocities of the parts resisted. The general equations 
of motion will then be of the form 

&c. = 
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To solve these equations we write 

Substituting we obtain a determinantal equation to find m. 
If we put 

A^jy + A^m + A,, 5 = By + B,m + B; 
A' = ^>» + Aim + a;, &c. 

this equation may be written in the simple form 



A, B, C, 
A\ B, C\ 



= 0. 



We may also write the equation in the form/(m) = 0. 

The coefficients Jf, Jf', &c. are not independent, but if we 
represent the minors of A, B, C, &c. by a, b, c, &c. we may easily 
show that 



We also have 



a b c 






It may be shown by properties of determinants that these 
equations all give the same ratios. If A, be the second minor 
obtained from the determinant /(m) by omitting the first and 
second rows and columns, we know that 

A,/ (m) = aV — a'6. 

Hence if /(m) =: we have -7 = 77 . 
•^ a b 

In the same way we may show that -7 = - , and so on. This 

Qf 

property of Determinants is given in Dr Salmon's Higher Algebra, 
Lesson iv. Ex. 1. 

The general solution of the equation may therefore be written 
in the form 

z = &c. J 

where ij, L^... are arbitrary constants, a^, a^, &c. the values of 
the minor a when m^, m, . . . are substituted for m ; 6^,6^... the 

1—2 
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values of the minor b when similar substitutions are made> and 
so on. 

4. We see that the whole character of the motion will de- 
pend on the signs of the quantities m^, m,... If any one be real 
and positive, x, y, &c. or some of them will ultimately become 
large, and the steady motion about which the system is oscillating 
will be unstable. If all the roots are real, negative or zero and 
unequal, the motion will be stable. 

If two of the roots be imaginary we have a pair of imaginary 
exponentials. If these imagmaiy roots be a ± )8V — 1, the terms 
can be rationalized into 

The motion will be stable if a be negative or zero, and unstable 
if a be positive. 

If two roots be equal, the form of the solution is changed. Let 
wi, = mj + A where h will be ultimately zero, we then have 

If we now make L. and Z. infinite in the usual manner, we 
find 

&c. = &c., 
where Jfj, M^ are two arbitrary constants which replace Z^, L^. 

In the same way if three roots are equal we have 

This rule will be found convenient in practice to supply the 
defect in the number of arbitrary constants produced by equal 
roots. At present we are only concerned with their effect on the 
stability of the system. The terms which contain ^ as a factor 
will at first increase with ^, but if m be negative, the term t^e^ 



can never be numerically greater than — . If m be very small 



n 
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the initial increase of the terms may make the values of x and y 
become large, and the motion cannot be regarded as a small 
oscillation. But if the system be not so much disturbed that 

M. — is large, the terms will ultimately disappear and the 

motion may be regarded as stable. If, however, the real parts of 
the equal roots are positive or zero, the terms will become large 
and the motion will be unstable. 

6. In some cases, however, the relations which exist between 
the coefficients are such that the terms which contain ^ as a factor 
are all zero. It is of some importance to discriminate these cases, 
for the stability of the system is then imaffected by the presence 
of equal roots. 

Let us suppose first that the determinantal equation has two 
roots only equal to m^ and let the terms depending on these be 

&c. = &c. 

Substituting in the equations of Art. (3) we have, following 
the same notation as before, 

^'iv;+5'iv;+c7Jv;"+... = o[ i. 

&c. = 3 



ji. 



&C. = &C, 

To avoid entering more minutely than is necessary into the 
properties of linear equations, we shall assume that these equations 
for the given value of m lead to but one solution with two of the 
^s arbitrary, unless the determinantal equatiou has more than 
two roots equal to m^. If in this unique solution the N^^ are all 
zero we must have 

^'ivr + /?iv;' + ... = of III. 

&c.=0 
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Since two of the constants N y N\ &c. are to be arbitrary, let 

. N N' 
them be N,. N' then since —^=-jL we must have the minors a 

N N'' 
and b each equal to zero. Also since -^ = — ^, we shall have N^' 

infinite unless c = 0. In the same way we may prove that all the 
other first minors are zero. 

And if the first minors are zero, we may show that two of the 
equations may be deduced from the others. Let the symbol 

J I TV represent the second minor, with the usual sign, formed 

)y omitting the rows and columns in which AB A'B' occur. 

Then since the minors a, 6, c, &c. are zero, we have 



^U4]+^''U'^'j+---=^ 



JV. 



&c. = 0. 
Omitting the first line of III. let us multiply the others by 



VAB^VAB 
[A!Bf\ \A!'ff'_ 



&c., respectively. Adding the results, we have 



an identity. Hence the second equation may be deduced from 
the others which follow it. In the same way, the first equation 
may be deduced from the others. 

Rejecting the first two equations, let us transpose the arbitrary 
constants N^ and N^ to the right-hand sides of the remaining 
equations. If there are to be only two arbitraiy constants, these 
remaining equations must be independent; solving, we have 

\_A'B'\ -"' ^^' [eG'\ + ^^ [A'CX 

with similar equations for the others. Hence the constants 
N^, N^, &c., are connected by equations of the form 

so that when any two are chosen as the arbitrary ones, the others 
may be deduced from them. 

If the determinantal equation has three roots equal to m^, and 
if the terms which contain < as a factor are all zero, the equa- 
tions III. must admit of a solution with three of the constants 
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N^, N^\ &c. arbitrary. If these be iV^, N^\ J^/', we see that the 
second minors jypf\ a'P'\\ A'jy\ ^^st be zero. 

Next since 

we must have d/t^J and jfrvA both zero or N^'^ infinite. 

Hence all the second minors are zero. And if the second 
minors are all zero, we have 



I// 



■ ABG 
A'JffC 
A"B'G" 



+ A 



/// 



■ ABG 
ABG' 
A"S"G 



tn 



+ ...F:0, 



and by similar reasoning three equations may be deduced from 
the remaining ones. We have then 



N. 



' BGD ' 


-N' 


' AGD ' 


+n:' 


r ABD 1 


-Nl" 


' ABG ' 


B'G'Df 




A'G'B 




A'BU 




A'B G' 


_B'G"iy'. 




A" a'B' 




rs'if'. 




A"B" G'. 



=0, 



with similar equations. 

Since f{rn) is the determinant formed by eliminating N^^ N^\ 
&c. from III. we have 

■ • 

dfjm) ^df(m)dA ^ df(m)dB ^ ^,^ dfim)dA' ^ ^^^ 
dm dA dm dB dm * dA! dm 

= a-j — h6j— + &C. + a' -5 — h &c. 
dm, dm dm 

This vanishes when a, J, &c., a', &c. are all zero. If therefore the 
first minors of /(m) all vanish when w = m^, the equation /(m) = 

has two roots equal to m^. In the same way -y- vanishes if all its 

first minors are zero. But 






drri 



dtti^ dm dm 



da 



vanishes if a, -v— , &c. are all zero. If therefore the first and 

dm 

second minors of /(m) all vanish when m = Wj, the equation 
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fipi) =0 has three roots equal to m^. It is evident the proposi- 
tion may be extended to any number of roots. 

The test that, when equal roots ocdur in the determinantal equa- 
Hon, the terms in ilie values ofx, y, &c. which contain t cw a factor 
should be absent may be stated thus. If tliere are two equal roots 
all the first minors must vanish. If three equal roots, all. the first 
and second minors must vanish, and so on. In these cases the 
equal roots introduce merely a corresponding indeterminateness 
into the coefficients. 

When there are more equal roots than there are rows in the 
determinantal equation, it is easy to see that there must be some 
terms in the integrals which contain ^ as a factor. 

[The following simple example will illustrate the application 
of this test 

A particle is in equilibrium at the origin of co-ordinates under 
the action of forces whose force function U is given by 

If the level surfaces ai-e ellipsoids and the force acts inwards, 
it is clear that the equilibrium of the particle must always be 
stable. If then any equal roots occur in the determinantal 
equation^ the test should show that the terms which contain t as 
a factor are absent 

If T be the semi vis viva of the particle and if its mass be 
taken as unity, we have 

Omitting accents and forming the discriminant of — mT+ Z7 
we have the following determinantal equation : 



A-m^ F E 

F B^m' D 
E D G-m^ 



= 0. 



This is the "discriminating cubic" which determines the axes 
of the quadric Z7= c, where c is a constant. The conditions that 
two of its roots should be equal, i.e. that the quadric should be 
a spheroid, are well known to be 

.. EF ^ FD ^ DE , 

where m^ is equal to either root These are just the conditions 
obtained by equating any first minor of the determinant to zero. 
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The conditions that three of the roots of the cubic should be 
equal, i.e. that the quadric should be a sphere, are 

These are the conditions that every second minor should vanish. 

In this example we have taken the case of a single particle. 
Similar remarks however apply when any system of bodies is 
disturbed from a state of stable equilibrium. The oscillations 
may be found by the method of Lagrange. The final determi- 
nantal equation may be conveniently formed by equating to zero 
the discriminant of — mT+ U, where Tis the semi vis viva with 
the accents denoting diflferentiations with regard to the time 
omitted, and U is the force function. It is a known theorem that 
the existence of finite equal roots does not affect the stability of 
the equilibrium. Hence the conditions for equal roots must be 
such as to make all the minors equal to zero. Conversely, this 
theorem will often conveniently give the conditions that Lagrange's 
determinant has equal roots.] 

6. That there should be a difference in the modes in which 
equal roots affect the motion is no more than we should expect 
a priori. Suppose the coeflBcients of the equation /(m) = to be 
functions of some quantity n, and that as n passes through the- 
value tIq, two roots become eaual to each other. Let the quadratic 
factor containing these roots be m' + 2am 4- /8, and let us consider 
only the case in which a and fi are real. We have a' — )8=0 
when n = nQ. If a* — /8 change sign as n passes through the value 
n^ the roots will change from a trigonometrical to a purely ex- 
ponential form, which would indicate a change from oscillatory to 
non-oscillatory motion. The passage from one kind of motion to 
the other may be effected through a motion represented by expres- 
sions having the time as a factor. But if a' — /8 does not change 
sign, for example, if it be a perfect square for all values of n, there 
will be no change from one kind of motion to the other, and in this 
case we should expect that the motion when the roote are equal 
will be represented by terms of the same character as before. 
Briefly, we may expect equal roots to introduce terms with ^ as a 
factor at the boundary between stability and instability ; and to 
introduce merely an indeterminateness into the coeflBcients when 
the motion is stable on both sides. 

It is easy to show that in the first of these two cases the 
minors could not contain either of the factors of m^ + 2am -I- )8. 
For since c? — P changes sign, these factors are in one case ima- 
ginary ; and therefore if one fector occur in any minor the other 
must also be present. The minors would not only vanish, but 
must have equal roots also. But as in Art. (3), \f{m) ^aV —ah. 
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Hence if all the first minors have equal roots it is clear that either 
/(m) has more than two equal roots, or all the second minors must 
vanish. The latter is impossible unless /(m) has more than two 
equal roots. 

These general considerations are not meant to replace the 
proofs given in the last article, but merely to explain how a 
difference in the effects of the equal roots might arise. 

7. Summing up what precedes, we see that if a dynamical 
system have n co-ordinates its stability depends on the nature of 
the roots of a certain equation of the 2?ith degree. 

If the roots of this equation are all unequal, the motion will 
be stable if the real roots and the real parts of the imaginary roots 
are all negative or zero, and unstable if any one is positive. If 
several roots are equal the motion will be stable if the real parts 
of those roots are negative and not very small, and unstable if 
they are negative and small, zero, or any positive quantity. But 
if, as often happens in dynamical problems, the terms which con- 
tain ^ as a factor are absent from the solution, the condition of 
stability is that the real roots and the real parts of the imaginary 
roots of the subsidiary equation should be negative or zero. 

8. When the equation /(2>) = is of low dimensions we may 
solve it or otherwise determine the nature of its roots ; the 
stability or instability of the system will then become known. 
But if the degree of the equation be considerable this is not a very 
easy problem. We shall devote the two next Chapters to the 
consideration of two methods by either of which, without solving 
the equation, we can determine the conditions that the real roots 
and the real parts of the imaginary roots should be all negative. 
The determination of these conditions has, it appears, never before 
been accomplished.* The consideration of the equations of motion 
will then be resumed, and the form of the determinantal equation 
/(Z>) = when the forces admit of a force function will be more 
particularly investigated. 

* [These conditions for the cases of a biquadratic and a quintic had been found 
by the Author in 1873, and read before the London Mathematical Society in June, 
I87I See also the third edition of the Author's Rigid Dynamics, Art. 436.] 
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Statement of the theory/ by which the necessary and sufficient tests of 
stability are found. Objections to this theory. Arts, 1 — 3. 

These tests shown to be integral functions of the coefficients. Art, 4. 

Method of finding these tests when the coefficients of the equ>ation a/re 
numerical, or wfien several terms are absent. Arts. 5, 6. 

All these tests sJiown to be derivable from one called tfis fundamMnUd 
term. Arts, 7, 8. 

The fundamiental term found as an eliminant. Art, 9. 

A method of finding the fundamental term by derivation from the funda- 
mental term of one degree lower. Art, 10. 

Another a/nd better method of doing the same by means of a differential 
equation. Arts, 11 — 13. 

1. The object of this Chapter has been explained at the end 
of Chapter I. Briefly, the criterion that the motion of a system 
of bodies should be stable is that the roots of a certain equation 
should have all their real parts negative. We propose to investi- 
gate these conditions. 

Let the equation to be considered be 

Let the real roots be a^, a,... and the imaginary roots be 

a,±AV-l, a,±/S,\rn, &c. 

Then 
f{x)=p, (a?«aj {x-a,) ... {a? - 2ol,x + a^ -¥ P^, &c. 

If then ttj, a., &c. a^, a^, &c. are all negative, every term in 
each factor, and tnerefore in the product, must be positive. 

It is therefore necessary that every term in the ecjuation 
/(a;) = should have the same sign. It will be convenient to 
suppose this sign to bo positive. 
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It is also clear on the same suppositions that none of the 
coe&cienis Pq, p^, >».p^ can be zero, except when the roots of the 

equation are all of the form i/SV— 1, or when some of the roots 
are zero. 

2. Let us now form the equation whose roots are the sums of 
the roots oif{x) taken two and two. Let this be 

n — 1 

where m = n — a— . The real roots of this equation will be 

aj + a,, tti + a,, &c. 20^, 2afj, &c. and the imaginary roots will be 

a^ + ttj + /8j V — 1, &c. It is clear from the same reasoning as 
before that if a^, ttj.&c. a , or,, &c. are all negative, the coefficients 
P^ Pj, &c. must all have the same sign. 

Conversely, if ^o> Pi • • • l^^ve all the same sign, the equation 

f(x) can have no real positive root, and if P^,, P^ ... P^ have all 

the same sign the equation F{x) can have no positive root, and 

therefore f{x) can have no imaginary root with its real part 

positive. 

3. Our first test of the stability of a dynamical system is that 
all the coefficients of me dynamiccd equaticm f (D) = and all the 
coefficients of its derived equation F (D) = should have the same 
sign. 

It should be noticed that though these conditions are all 
necessary and sufficient, they are not all independent. We obtain 
too many conditions. In many cases, however, we can at once 
reduce them to the proper number of independent conditions, and 
when this is difficult we can have recourse to the second method, 
to be given in the next Chapter, which is free from this objection. 

In order to apply this method with success, it is necessary to 
have some convenient methods of calculating the coefficients 
P P P 

4. The first method which suggests itself is one similar to that 
usually given to determine the coefficients of the equation whose 
roots are the squares of the differences of the roots of any given 
equation. 

. If ^1, 5j... be the sums of the first, second powers, &c. of the 
roots of the equation /a? = 0, we have by Newton s theorem 

where p^ has been put equal to unity. If 2^, 2,... be the sums of 
the powers of the equation F{x) = 0, we have in the same way 

2 4-P2 +P2 + =0- 
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we may also prove 

2,= («-l)S., 

S.= (»-4)5.+ 3^.^„ 

2. = (n- 8) Si + 45,5, + 3/8'/, 

2, = (n- 16) 5. + 55^ + 105,5.; 

and the general relation can be found without difficulty. 
In this way we find 

p_ («-l)(n-2)(n-3)(«-4) , («-2)'(n-3) 
♦ 1.2.3.4 P^^ 1.2 P'^* 

+ (n-3)>.y. + ^""f/^~^^ J)/+(n-8)y.. 

But the process becomes longer and longer at every stage. 
We shall therefore proceed to point out some other methods of 
obtaining the coefficients. 

This method of proceeding has indeed been stated only because 
it proves in a convenient way that whei^ p^ = 1, all the coefficients 
P^f P^, P^... of the derived equation are integral rational func- 
tions of the coefficients p^, p^. . . p^. 

6. The equation f (x) = beinff given, to calculate the coefficients 
o/F(x)=0. 

Put x^y ±z and equate separately to zero the sums of the 
even and odd powers of z, we have 






/(y)+/"(y)^+/"(y)g+..- = o 



a» 



> • 



/'(y)«+/"'(y)g + ... =0J 

Rejecting the root -? = 0, let us eliminate z. Then the roots 
of the resulting equation in y are the arithmetic means of the 
roots of/(a?) = 0. 
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If, on the other hand, we eliminate y we have an equation of 
an even degree to find z. This, putting 4«'=f, is the equation 
whose roots are the squares of the differences of the roots of the 
given equation. 

It may be thought that this elimination may prove tedious, 
but it will be presently shown that only the first and last terms of 
the result are really wanted. All the others may be omitted in 
the process of elimination, and thus the labour will be greatly 
lessened. The method is however most useful when the given 
equation has several of its terms absent. 

6. Eocample. To determine the condition that the roots of the 
biquadratic 

X* +px' + qx* + rx + s = 
should indicate a stable motion, . 

Applying the rule we have 

+ (2p^q+pr + g*— 45) 0^"+ (pg* +^V- 4^5) oj+pgrr — V-^V 

The first four coefficients contain only positive terms, and need 
not be considered. If the last three coefficients be called P4, -Pg, 
P^, we have 

pP,-^P,=^{pq-2ry + p'r. 

If then Pg is positive, all the other coefficients are positive. 

The necessary and sufficient conditions of stability are there- 
fore that p, q, r, s should be finite and positive, and 

Pg =pqr — r' —p^s 
positive or zero. 

7. In forming the derived equation F(x) the only difficulty 
is to form the last term P^. For when this is known the other 
terms can he at on^ce derived from it by an easy process. 

Let a, b, c... be the roots of f{x) = with their signs changed, 
and let 

Let A stand for the operation 

da db dc 
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Then since — = abc ... we have obviously 
In the same way we have 

^ Pn-l _ O Pn-9 

Po Po ' 



And generally 



and so on up to 



Po Po 



Po 



Let us now operate with A on any expression 

which has the same number of factors in every term. Let r be 
the number of factors, then <f> may be written 

= {m ^+(n- l)p. ^^ + ... +p^, 1^} f 

Let Pq, P^...P^ be the coefficients of the derived equation 
F(x), and let F^ = 1. Then since the roots of F(x) are the sums 
of the roots oif{x) taken two and two, it is easy to see that 

3P_. = iAP_^ 

&c. = &c. 

Thus when P^ is known, the other terms may be calculated 
without difficulty. The term P^ will be called the fundamental 
term of the equation. 

EocampU. Given in the case of a biquadratic 



to calculate P^. 
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Performing the operation 
on Pj we find after division by 2, 

-P5 = PiP* + PiPs " ^,PiPv 

which is the result already given. 

8. It should be noticed that in the equation 

/(^) =|)X +l>ia:*-^ + ... +i>n = 0, 

if we regard ^ as a number, p^^p^ ... 2>n ^-re all of equal dimensions. 
It follows from the theory of dimensions, that if any subject of 
operation be the sum of a number of terms of the form 

there must be the same^ number of factors in every term. For 
example, in every term of the expression for P« we have 
a + iS + 7 &c. the same. 

On the other hand, we may regard a? as a quantity of one 
dimension, and in this case jPo>J^i "*Pn ^ave their dimensions in- 
dicated by their suffixes. We must therefore have /8+ 27+88+... 
as well as a + /8 + 7+ ... the same in every term. 

These two tests of the correctness of our processes will be 
found convenient. 

9. 7^e whole derived equation being knoton when the funda- 
rH'Cntal term is known, it is required to find the fundamental term. 

[First Method.] 

If we write — a? for a? in any equation, we have a second equation 
whose roots are equal and opposite to those of the first equation. 
If we eliminate x between these two, we shall get a result which 
must be zero when the two equations have a common root. The 
eliminant must therefore contain as a factor the product of the 
sums of the roots of the given equation taken two and two. 

It will afterwards be shown that the last term P^ of the 
derived equation (when p^ is put equal to unity) always contains 
the term 

PvPi^Ps-' 
with a coefficient which is positive and equal to unity. 
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Hence we have this rule, to find P . eliminate x between 

or +i)X"^+jpX"* + ... = 0' 

and divide the result by the coefficient o{pj)j>^,,.p^^^. 

It is obvious that the result may be written down as a 
determinant. On trial, however, it will be found more convenient 
to make the elimination by the method of eliminating the highest 
and lowest terms than to expand the determinant. 

10. Given the fundamental term of the equation derived from 
/{x)=p^''-'+p,x''--+... + p,., = (1), 

to find the fundamental term of the equation derived from 

i>o^*+;>xa?*-*+...+l>H-i»+i)» = (2). 

[Second Method.] 

Let Q^ be the product of the sums, two and two, of the roots 
of the equation (2) taken with their signs changed, so that Q^ is 
the same as the fundamental term of the derived equation and 
diflfers from P^ only in having a suffix more convenient for our 
present purpose. 

Let Q^ be expanded in a series of powers ot p^ : thus 

where (l>^ (f>^, ^j, &c., are all functions o{p^ »j, &c., which functions 
have to be found. Let the roots of /(a?) = with their signs 
changed be a, b, c..., then 

^».i=(a+5)(a + c)... 

Let us introduce a new root, which, when its sign is changed, 
we shall call r. Then 

Q^={a + b){a + c) ... (r + o)(r + i)... 

This value of Q^ must be the same as that given by the series 
when we write 

Po> Pi + ^Po* i>s + »yi» &c., l>,.i + rp..,, rp„.^, 

respectively, for p^, p , p^, &c., p„_^, p^. Equating the coefficients 
of the terms independent of r we have 

R. A. 2 
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Equating the terms containing the first power of r we have 

Substituting for ^^ we have, by a known theorem in the Diflferential 
Calculus, 

Equating the terms containing the second power of r we have 






and so on. Thus we have 

If we examine this process, we see that when Q„__^ is known, 
we may at once write down the terms independent of p^ and the 
coeflBcient of p^. The process to find the coeflScient of p^ is 
longer, but it may be much shortened by the consideration that 
when Po = 1 the result must be an integral function of the coeffi- 
cients. We may therefore omit all terms as soon as they make 
their appearance, which do not contain the factor |>Vi> ^^^ we 
know that such terms must disappear from the result. 

This method is not so convenient as the one which will be 
presently given to find the coefficients of the higher powers of p^. 
But it is useful as showing that Q^ contains the term 

with a positive integral coefficient equal to unity. This will be 
clear from the consideration that the term independent of p„ in 

Q^ is obtained from Q^.^ by multiplying by ^-^^ . If therefore Q^^ 
contains ^'^''j^lf""' , Q^ must contain the term ■ P'^g'^;^*^^ . No 

Po , , ^ ^ Po 

other term can be formed which is equal to this with an opposite 
sign, for the terms which enter by the other processes to be per- 



ir.] THE THIRD METHOD. 19 

formed on Q„_^ all contain |>» as a factor. Now Oj= - ; therefore 
Q^ contains the tenn ^-^ , ft contains ^^'g , ^^d so on. 

Po Po 

It has been shown that every term in p^'^Q^ has the same 
number of factors (Art. 8). It follows from this reasoning that 
this number of factors is w — 1. 

11. To find the fundamental term of the derived equation hy 
means of a differential equation, 

[Third Method.] 
The fundamental term required is a factor of the eliminant of 

Let a;* = y, then we have 

py'+p»y*"+ =o) 

i'.r+P.y'^+ =0}nodd. 

n-1 n— 3 I 

PiV' +i>8y" + = ^' 

If we write y + dy for y the result of the elimination must be 
the same. Hence if we make 






even, 



n odd, 



dy 2 P"' dy- 2 ^^'^■ 

rfy 2 ^" dy 1 ^•'**'J 

and if ^be the eliminant, we have 

dy 

It follows that whether n be even or odd, J? must satisfy the 
equation 
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dE dE 

„ / dE dE\ 

, _ / dE dE\ 

+ 

We may make the elimination by multiplying the two equa- 
tions by y, y'..., until we have as many equations as we have 
powers to eliminate. If in the determinant thus formed, we 
multiply out the terms in the diagonal joining the right-hand top 

——1 — +1 
corner to the left bottom comer, we get when n is even p " « * 

and when n is odd p^^ p^^ - Now Q^ must contain n — 1 factors 

w — 1 

and be of the n — ^ — th degree. Hence when n is even j&= opji^ 

and when n is odd j&= cQ^ where c is some constant. 

Now -T— does not occur in the above diflferential equation. 

Hence treating p^ as a constant, we see that Q^ must satisfy the 
differential equation 

V^"- rfp. ^"- dp J + l^- dp^, + P-^ dp J 



+ 3(p^f^ + p^4^+&c. = 0. 



12. We may show that pj^'^ Q^ is a symmetrical function of the 
coefficients p^, Pj..., Pn and the same coefficients read backwards. 
Let a, 6, c... be the roots of f(x)=0 with their signs changed, 
then <?n = (^ + ^) (a + c).... 

If now we read the coefficients in the opposite order, the roots 
of the equation thus formed will, when their signs are changed, be 

-,=-.... If 0'^ be the fundamental term of the equation derived 

a o '• 

from this, we have 



'^Hhl){l^\)-- 



Since — = ahc... we see that 

Po 

rr^Qn^Pf^^'^Qn' 
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We may therefore infer that ^^*~*Q» also satisfies the diflferential 
equation 

13. We may use either of these differential equations to find 
Q» when Q^^^ is given. 

Let the first diflferential equation be represented by 

and let 

where A^y A^y ... are functions oi p^, p^, &c. The value of A^ has 
been proved in Art. 10 to be 

To find the other coeflScients of the powers of p^^ substitute 
this value of Q^ in the diflferential equation ; we have 

= VA^ +PnVA -^Pn-^^i 

+ &C. 

Equating the several powers of ^„ to zero, ve find 

&c. = &c. 

Thus by one regular and easy process each term may be 
derived from the other. 

In performing this process we may omit every term in the 
subject of operation which does not contain /?„_j. For j9„_j can be 

introduced only by performing -y— , and since p» is absent from the 

P* V . 

coeflScients, this operation yields nothing. 
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lu this way we find 

PoQs=PJ>»-PoP>. 

PoQ, = PlP»P, -PoP» -Plp4. 

Po Qt =PtP»P»P4 - P»P,'P4 - PxP! 

-Pt (-PoPiP, +PiP»-^PoPJ>^ 

+ 3^ (- 2a')- 

To illustrate this process, consider how Q^ is obtained from Q^, 
The first line is formed by multiplying the line above by p^, this 
is A^. To find the coefficient of —p^ we operate with 

on such of the terms in the line above as contain p^ and then 
divide by p^» Performing the same operation on the coefficient of 

("Pb) we obviously obtain the coefficient of f-^ . 

In M. Serret's Cours cFAlghbre Sup^rieure,'Note III., there will 
be found a method of forming the last term ojf the equation to the 
squares of the diflferences, which suggested the method used in 
Art. 13, of substituting in a differential equation, if only a differ- 
ential equation could be found. [See also Dr Salmon's Higher 
Algebra, Arts. 60, 64?, and 72.] ' 



CHAPTER III. 



Statement of the tlieort/ by which the conditions of atahitity of a dynamical 
system with n co-ordinates are made to depend on 2a conditions. 
Arts, 1 — 4, 

A rule by which these conditions may be derived one from anoth^r^ to- 
gether with certain other true but not independent conditions. Arts, 

A rule by which, when the coefficients of the dynamical equation a/re letters, 
the 2n conditions of stability m^ay be inferred, one from another, 
without writing down any other conditions. Arts, 9, 10. 

A mstliod by which certain extraneous factors m>ay be discovered and 
omitted. Arts, 11, 12. 

Consideration of the reserved case in which the dynamical equation has 
equal and opposite roots. Arts, \^ — 18. 

A few geometrical illustrations not necessary to the argument. Arts, 
19—25. 

Application to a Dynamical Problem, Art, 26. 

1. It has been shown in the first Chapter that the stability of 
a dynamical system with n co-ordinates oscillating about a state of 
steady motion depends on the nature of the roots of a certain 
equation of the 2n^ degree which we may call 

The system is stable if the real roots and the real parts of the 
imaginary roots are all negative. Now Cauchy has given the fol- 
lowing theorem of which we shall make some use. 



Let « = a? + y V— 1 be any root, and let us regard x and y as 
co-ordinates of a point referred to rectangular axes. Substitute 
for z and let 

Let any point whose co-ordinates are such that P and Q both 
vanish be called a radical point. Describe any contour, and let 
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a point move round this contour in the positive direction and 

P 

notice how often ^ passes through the value zero and changes its 

sign. Suppose it changes a times from 4 to — and 13 times from 
— to +. Then Cauchy asserts that the number of radical points 
within the contour is ^ (a — /3). It is however necessary that no 
radical point should lie on the contour. 

2. Let us choose as our contour the infinite semicircle which 
bounds space on the positive side of the axis of y. Let us first 
travel from y = — oo toy = +oo along the circumference. 

we have changing to polar co-ordinates 

f{z) =py (cos nO + smnO ^/ — I) + , . , 

Hence 

P=j9or"cosn^ + i>ir""*cos («-l) ^+... ) 
Q = p^r"" sin nd + p.r*^* sin (n-^ 1) ^ + ... j' 

In the limit, since r is infinite, 

j^ = cot nO ; 

^ vanishes when nO = (2/c + 1) —, i.e. 

0=±-'l, ±-^ ±-? (A); 

7i2n2w2 ^^' 

pr is infinite when n& = 2k ^ , i.e. 

*-«. 41 4h 4i («)• 

The values of 6 in series (B) it will be noticed separate those 
in series (A). 

P . 

When is small and very little greater than zero, ^ is posi- 

tive, and therefore changes sign from + to — at every one of the 
values of in series (A). If n be even there will be n changes 
of sign. If n be odd there will be n — 1 changes excluding 

0= ±^ i in this case ^ is positive when 5 is a little less than — , 
and negative when ^ is a little greater than j^ . 



and 
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Let us now travel along the axis of y still in the positive di- 
rection, viz. from y=4-oo toy=--oo. Since ^ = it will be 
more convenient to use Cartesian co-ordinates, we have, since 

and « = yV— 1, 

The condition that there should be no radical point within the 
contour is that this expression should change sign through zero 
from — to + as often as it before changed sign from + to — on 
travelling round the semicircle. If n be even the numerator has 
one more term than the denominator, and when p^ and p^ have 

the same sign, J begins when y is very great by being negative. 

In order that it should change sign through zero n times, it is 
necessary and sufficient that both the equations 

Pn -Pn^f +Pn.y - • • • = 0, 
P»^iy - Pn-s^* +Pn-6^ - . . . = 0, 

should have their roots real, and that the roots of the latter should 
separate the roots of the former. 

If n be odd, the numerator and denominator have the same 

P 
number of terms, and when p^ and p^ have the same sign, ^ begins 

when y is very great by being positive. In order that iT should 
change sign through zero from — to + n — 1 times, it is necessary 
and sufficient that the same two equations as before should 
have their roots real, and that the roots of the former should 
separate the roots of the latter. 

In order then to express the necessary and sufficient conditions, 
that f{z) = may have no radical point on the positive side of 

the axis of y, put z = y V— 1 and equate to zero separately the 
real and imaginary paHs. Of the two equations thus formed, the 
roots of the one of lower dimensions must se^mrate the roots of the 
other,' It is also necessary that the coefficients of the two highest 
powers ofz in f (z) should have the same sign, 

3. It has been stated that p^ and p^ the coefficients of the two 
highest powers in f{z) must have the same sign. It is easy to see 



26 EXTENSION OF STURM'S THEOREM. [CHAP. 

P 

that, if they had opposite signs, ^^ would change sign through zero 

2/1 times as we travel round the contour. All the radical points 
Df the equation would then lie on the positive instead of the 
negative side of the axis of y. 

It has also been assumed that no radical point lies on the 
contour. It has therefore been assumed that/(;5) = has no root 
of the form « = ly v — 1. It will be more convenient to consider 
this exception a little further on. 

4. It is required to express in an analytical form the conditions 
that the roots of an equation f^ (x) = may be all real, and may 
separate the roots of another equation f^ (x) = q/* one degree higher 
dimensions. 

To eflfect this, let us use Sturm's theorem reversed. Perform 
the process of finding the greatest common measure of f (x) and 
f^ {x)y changing the sign of each remainder as it is obtained. Let 
the series of modified remainders thus obtained be /g(a;), y^ (a;), 
&c. Then it may be shown that when any one of these functions 
vanishes, the two on each side have opposite signs. It is also clear 
that no two successive functions can vanish unless /, {x) and f^ [x) 
have a common factor. This exception will be considered pre- 
sently. 

Hence in passing from a? = — oo to + oo no variation of sign 
can be lost except when f (x) vanishes. If a variation is lost it is 
regained when x has tbe next greatest value which makes f (x) 
vanish unless :f/(x) = has a root between these two successive 
roots otf^ix) =='0;; Hence this rule: — 

The roots of .the equations fi(x) = 0, ^(x) = 0, will he all real 
and the^rpots of the latter will separate those of the former, if in 
the series 

fi(x), faCx), faCx)... 

OS many variations of sign are lost in passing from x = — oo to 
X = H- 00 as there are vmts in the degree of the equation f^ (x) = 0. 

We have supposed the variations of sign to be lost instead of 
gained in passing from ic = — oo to + oo . That this may be the 
case the signs of the highest powers of f (x) and f^ (x) must be 
the same. 

These functions are alternately of an even and odd degree, the 
condition that the whole number of variations of sign may be lost 
in passing from a? = — oo to + oo may be more conveniently ex- 
pressed thus : — The coefficients of the highest powers of x in the 
series 

fiW, fa(x), f,(x) ... 

mv^t all have the same sign. 
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5. The proceaa of finding the greatest common measure of 
two algebraic expressions is usually rather long. We may in our 
case shortea it materially by omitting the quotients and perform- 
ing the division in the following manner. Let 

/i i'^) =P<P'' -Pj^'' + pX"* - ■ ■ ■ 
/. ('^) =Pi^'" -P.^"* +J>.*"^ - ■ ■ ■ 
then, since p^ is positive, it easily follows by division that 

- f,{x) =Ax'^ ~ A'a^-^ -^ A"x°-^ - ... 
where A =PiPa—pj>ai 

^'^P.Pt-PoP,* 
&c. = &c., 

80 that by remembering this simple croasrmultiplication we may 
write down the value of f, (x) without any other process than what 
maybe performed by simple inspection. In the same way^(j;), 
&c. may all be written down. 

6. Ex. 1. Express the conditions that the real roots and 
real parts of the imaginary roots of the cubic 

may be all negative. 

ft {^ = ^ — 5-c, 

The necessary conditions are that 

p, p2 — r and r 
must all be positive. 




Ex. 2. 
quadratic 



Express the corresponding conditions for the bi- 

/■W-K -'■■»'. 

/.W-Cm-')*" -p». 
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The conditions are that 

p, pq — r, {pq — r) r—p^s and s 
must be all positive. 

These are evidently equivalent to the five conditions that 

P» ?, ^, s, {pq-r)r-p\ 
should be all positive. 

In both these examples all the numerical work has been 
exhibited. 

7. Since the coeflScients of the highest powers of x in j^ (^) and 
f^{oc) are p^ and p we see that the condition that p^ and p^ should 
have the same sign is included in the general statement that all 
the coefficients of the highest powers should have the same sign. 
If the function /(a?) be of n dimensions we thus obtain n necessary 
and sufficient conditions. 

On examining these conditions in the cases of the cubic and 
biquadratic it will be seen that they cannot be satisfied if any one 
of the coefficients of the given equation should be negative. 

8. Although the theorem in its present form gives n conditions 
as the proper number for an equation of the n^ degree, yet it is 
important to notice that it gives other conditions also which are 
true and may be useful. It has been shown in the second Chapter 
that all the coefficients of the equation f(x) = must be positive, 
hence the roots oif^{x) = must all be positive. It may be shown 
also, that the roots of each of the functions fAx)y f^{x), &c. are 
separated by the roots of the function next below it in order. 
Hence the roots of all these functions must be positive, and there- 
fore in every one of the functions the coefficients of all the powers 
must be alternately positive and negative and not one can 
vanish. If however fj^x) and fj^x) have one or more common 
factors some of the functions /3(a?),y^(a;), &c. will wholly vanish. 

9. When the degree of the equation is very considerable 
there is some labour in the application of the rule given in Art. 5. 
The objection is that we only want the terms in the first column, 
and to obtain these we have to write down all the other columns. 
We shall now investigate a method of obtaining each term in the 

first column from the one above it without the necessity of writing 
down any expression except the one required. 

We notice that each function is obtained from the one above it 
by the same process. Now 



\ 
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The first and second lines will be changed into the second and 
third lines by writing for 

Po» Vv P^ Pb^ ^^' 

the values 

Pv P1P2 -PoPs^ P^> P1P4 - PoPip &c- 

If then in any term of any function we make these changes, we 
obtain the corresponding term of the function next in order. 

10. Example. Express the conditions of stability for the 
quintic 

f{x) =py +;),«* +^X -hp^x"" + p^x +p,. 

We have 
/,(a:)=p^aj*+... 

f,(^) = {(j>lP2-PoPs)Ps-'Pl(PlP4-PoP5)}^> 

Ai^) = l[iPlP2 - PoPs)P3 - Pi (PiP4 - PoP^} (P1P4 - PoPd 

-{PiP^-PoPsY-P^^* 
fM = [{{PiP^ - PoP^P^ -P, {PyP^ - PM (P1P4 - PoP^ 

- iPlP2 - PoP^^ .i>J (PJ>^ -PoPs)P5' 

11. On examining the conditions as given in the cases of a 
biquadratic and quintic it will be apparent that several contain 
the previous conditions as factors. Thus the analytical expressions 
are rendered much longer than is necessary. It is now proposed 
to investigate a method of discovering and omitting these extraneous 

factors as they occur, and thus obtaining the required conditions in 
their simplest forms. 

Let the coefficients of the several powers of x in the functions 
be when taken positively 



fli«>)='p»> 


Pv 


i'4, 


P«--- 


fM = Pv 


Ps. 


Ptf 


p,... 


A{^)=A 


^', 


A", 


A"' 


/M = B, 


B', 


-B". 





&c. = &c. 
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Let us first find which of these terms contain j?i as a factor. 
Putting ^1 = and using the rule in Art. 5, the series become 



Pat 


Pv 


P4. 


p»' 


0, 


p»' 


P^' 


P,y 


PoP^t 


- PoPs> 


-PoP7> 


—PoP»> 


PoPi> 


-PoPiPt' 


-P<,P>P7> 


-PoPzP»> 


0, 


0. 


0, 


0. 


0, 


0. 


0, 

&c. 


0, 



Hence the C*s and D's all vanish and therefore contain ^^ as 
a factor. By the rule in Art. 5, the ^'s contain p^, the jP's con- 
tain p^, the G'sp^y the ITspi, and so on. 

But since each line is formed from the preceding by a uniform 
rule, it follows that the D's and ^'s contain -4 as a factor, the 
F*8 contain A\ the 0*s contain A^, the H's contain A\ and so on. 

The factor A in the Z)'s and E*s has its origin in the factor p^ 
which occurs in the (7s and B's and would not appear if that 
factor had been omitted when the G*s and D's were formed. The 
factor A in the D's and IPs in the same way gives rise to the 
factor B in the E's and F's, So that if we take care each time 
we perform the process described in Art. 9 to omit the common 
factor p whenever it occurs, all these subsequent factors will never 
make their appearance. 

We shall now show that if these factors are omitted, the 

w — 1 

dimensions of the n*^ function^ (a?) will be n — ^ — . First consider 

the actual dimensions of each function before the factors arc 
omitted. If we examine the rule by which each function is 
derived from the preceding, it will become evident that, the 
dimensions of each letter being indicated by its suffix, the 
dimensions of any function are equal to the sum of the two pre- 
ceding -1-2. 

In the following table the first column indicates the function. 
In the second column will be found the dimensions of the leading 
coefficient of that function when calculated by the rule in Art. 5. 
In the third column will be found the dimensions as given by the 

n — 1 

formula n — ^ . In the remaining columns are the dimensions 

of the extraneous factors ^j. A, B, &c. introduced into each term. 
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fM 














fM 


1 


1 








/sH 


3 


3 








/«(^) 


6 


6 




- 




/.(^) 


11 


10 


1 






/,(^) 


19 


15 


1 3 






XH 


32 


21 


2 3 6 






f.i^) 


53 


28 


3 6 6 


10 




/.W 


87 


36 


5 9 12 


10 


15 


&c. 


&c. 


&c. 


&c. &c. &c. 


&c. 


&c. 



&c. 

Each term in the second column is the sum of the two terms 
just above it + 2. The p!^ term in the third column is equal to 
the term just above it + (ri - 1). In all the other columns each 
term is the sum of the two terms just above it. The last term in 
the w*^ row is equal to the (n— 3)*** term in the third column. 
We wish to. show that any term in the second column is equal to 
the sum of all the terms in the same row to the right of that term. 
It is not difficult to show from the data just given that if this be 
true for any two adjacent rows, it is true for all the others, and 
hence we may assume it to be always true. 

It is clear that these extraneous factors may be omitted since 

by the conditions already expressed they are all positive. When 

omitted as they occur, the dimensions of the n^^ function has just 

n — 1 
been shown to be n — ^— . It is easy to see that the conditions 

thus reduced must contain the terms 

Pi* PiP,» PiP^Ps* PiP^PsPi^ &C. 

Now if we take any one of these as 

PiP^PsPiP^Pfi Wy 

and operate by the rule in Art. 9, we have 

(PiP^ - PoPs)Ps (PJ>i - PoP5)P, (PiP, - PoP,)pv 
which contains the term 



PiP2PsPiP4P5PiP6P, (2). 

Thus we have p^ introduced as often as there is a factor in (1) 
with an odd suffix. But it should be introduced only once. These 
extra p^a are the extraneous factors to be omitted. Each of these, 
if left, would appear as the fsictor p^p^— p^p^ in the next condition, 
and be still more complicated in the next after that. 
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In order then to obtain the several conditions in their simplest 
form it is only necessary after performing the operation described in 
Art, b or Art, Q to divide by /?j*, where k is one less than the number 
of factors with odd suffixes in the condition operated on, 

12. Example. Express the conditions of stability for the 
sextic 

f(x) =p^x' + p^a^ i-p^x' + p^x"" + p^a? + p,x + p,= 0. 
We have 

f i^) = (PiPiPs - PoPs - Pi Pi +PoPiP^) aj' + . . . 

A (^) = {PiP^PsPi - PoPsPi -PiP^ + ^P^PiPiPf^ -PiP^Pr^ +PoP2PzPf^ 

-P^P^ +Pi^P2Pe -PoPiPsPe] aJ* + . . . 

/e (^) = [PiP2PsPiP6 -PoPsP4P5 -PiP!p^+ "^P^PiPiP^-PiP^P^ 

+ PoPJ>sP6 - PoP' + ^PiP2P^P6 - SPoPiPsP^P, -PiP,PsPe 

■^PoPsPe+PiP^PiPe -PiPs] ^ + &C. 

/? (^) ^ coefficient of a? in ^ a? x by p^. 

13. In the preceding theory two reservations have been. made. 

1. In applying Cauchy's theorem it has been assumed that 
there were no radical points on the -axis of y. 

2. It has been assumed that P and Q have no common 
factor, so that none of the functions /j, /y, &c. vanish absolutely. 

If any radical point lie on the axis of y, it is clear that f{z) = 
must have a factor of the form (z^ 4- a^y. Let f{z) = («* + a*)*" <f>{z). 
In this case when we put ;2? =y V — 1, we have 

f(z) = {a'-yY(F+Q^/^y, 

.-. P=(a»-2/»)'P') 

Thus P and ^ have a common factor, and we are warned of the 
possible existence of radical points on the contour by the total' 
vanishing of some one of the functions /j,j^,/g, &c. 

The two reserved cases may therefore be included in one. If 
f{z) = be the equation furnished by dynamical considerations, 
we form the functions f^ f, f, &c. If all these be finite, the 
question of the stability of the system has been answered. If any 
one vanish absolutely, f (x) and f (x) have a common measure^ 
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and we must add some further considerations. It will be con- 
renient to examine separately the dynamical eflfect of the roots 
which do not and which do enter through the greatest common 
measure. Let us begin with the former. 

14. Following the same notation as before, we have 

Q = ±/a (y) =i>«-iy -i>«-8y' + • • •. 



If then/(«) have two roots, viz. ± (A + ft V — 1), which are equal 
and opposite, then /^ (y) and /, (t/) must have two common roots, 

viz. + — . The common measure therefore of f^ (y) and 

f^ {y) contains all the roots of f(y V — 1) which are equal and 
opposite. Conversely the greatest common measure of P and Q is 
necessarily an «ven function of y, and if it be equated to zero, its 
roots are necessarily equal and opposite. These roots must also 

satisfy/(yV"^) = 0. 

Let this greatest common measure be '^ («/') = 9, and let y*^ be 
the highest power which enters into it. Also let 

then <f} (f ) is a function which, as has just been shown, has not got 
two roots equal and opposite, and to this function we may apply 
Cauchy's theorem without fear of failure. Putting « = y V — 1, let 

P' 

Then we wish to express the condition that ^ should change 

sign from — to + through zero n — 2r times if n be even and 
w — 2r — 1 times if w be odd. But 



and /(y ^^=l) = ^{f){P' + q' V^); 

F 

Thus the number of changes of sign in jj is exactly the same 

p 

as that of -^ . The factor '^ (y') will run through all the functions 

/« (y)> f^ {y)» &c- obtained from /j (y) by a process which is equiva- 
lent to that of finding the greatest common measure of /j (y) and 

R. A. 3 
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/g {y). The changes of sign of this factor will therefore not affect 
the number of variations of sign in the series f^, /,, ^, &c. 

The last factor which is not zero isf^^_^^(]/) if w be the dimen- 
sions of/j (y). 

S.ence if we omit the considerations of the vanishing factors and 
a/pply the same rule as before to <A« n + 1 — 2r remaining factors, we 
can express the condition that the proper number of changes of sign 
from — to + have been lost through zero in the function (f> (z),i.e. that 
the roots not given by the vanishing of f^+i^ar ^^^ ^'^ ^f ^^ character 
to ensure stability, 

15. Let us next consider the effect on stability of the roots 
indicated by the absolute vanishing of one of the subsidiary func- 
tions. This function must be of the form 

where n is even. The corresponding factor of f{z) is 

^(^) = ?« + if^y + Jn-^^* + — • 

The roots of this equation are two and two equal with opposite 
signs, it is therefore necessary for stability that no root should 
have any real part. To express this condition, draw a straight 
line parallel to the axis of y at an indefinitely short distance from 
it, viz. x = h. Let us apply, in the same manner as before, Cauchy's 
theorem to the contour formed by this straight line and the infinite 
semicircle on its positive side. Putting z = h + z\vfe have 

F {z) =q, + 2q^Jiz' + q^/' + 4y„>/« + . . . 

The two functions are therefore, omitting the positive factor h, 

-P = ?n - ?iMi3^' + ?n-y - &c., 

Now/j (y) and/g (y) are what P and Q become when arranged 
in descending powers of y and the coefficients of their highest 
powers m^de to have the same sign. Hence 

-^'^^^ dy • 

The rule described in Art. 4 will now become the same as 
that usually called Sturm's theorem. We are to seek the greatest 
common measure of yj (y) and its differential coefficient, and make 
the coefficients of the highest powers of y in these two and in the 
series of modified remainders all positive. 



III.] OF A SUBSIDIARY FUNCTION. 35 

That we should have been led to Sturm's theorem in this case 
is just what we might have expected. For to express the condi- 
tions that the roots of 

are all of the form ± r V — 1 is the same thing as to express the 
conditions that the roots of 

are all real. 

16. There is however another mode of proceeding. Suppose 
we have calculated the functions /j, /,, &c. for the general equa- 
tion 

and find when the values of p^,, p^, &c. are substituted that some 
one function say f^ of the series absolutely vanishes, and therefore 
also all the functions which follow it. Then operate on each of 
these vanishing functions with 

repeating the operation until we obtain a result which is not zero. 
If we now replace these vanishing functions by these results we 
may apply the rule of Art. 4, just as if these were the functions 
supplied by the process of the greatest common measure. As this 
process is not so convenient as that already given it is unnecessary 
to consider it in detail*. 

17. As a numerical example, let us examine whether the roots 
of 

f(x) = a:*+ 2x' + 4aj' -j- 4a;' + 6aj*4- 6x» + 7aJ* + 4x -f- 2 =0 

satisfy the conditions of stability. In order to show the working 
of the method it will be necessary to exhibit all the numerical 
calculations. We have by Art. 5, 

• The fanction /^{x) vanishes because the eqaation /(«) = has two roots eqnal 
and opposite. If we pat z^gf + ht where h is as small as we please, this peculiarity 
will disappear. Thus if the valaes of z are of the form ^{a^p V-l) the corre- 
sponding values of z' are -hJ^a^p ^f-^' These values of z will indicate stability 
if a be zero and instability if a have any value positive or negative. If ^ be as smaU 
as we please and positive, the values of z' will indicate stability or instability under 
the same circumstances. We may therefore apply the rule of Art. 4 to the func- 
tion f{i^+h) instead of /(z), provided we retain only the lowest powers of h which 
occur. Hence all the functions /^(«), /j(«)-/r-i(«) which do not vanish are 
unaltered. To find what function will replace /,. (x) we must increase by ^ all the 
roots of /(z)=0 when their signs have been changed. This may be effected by 
performing on / (2) the operation represented b^ A in Art. 7. The rule in the text 
therefore follows from the one given in that article. 

3—2 
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/, {x) = 2aj' - 4a?' + Coj' - 4a?, 
/, (a;) = 4a?* - 6a:* + 1 Oa;* - 4, 
/4 (a?) = 4a:P - 4ai' + 8a?, 
/, (a?) = 8a?* - Sa;' + 16. 

Here we find/^ (a?) to be absolutely zero, accordingly by Art. 15 
we replace it by the diflferential coefficient oif^{x), this being 
3turm*s rule. We have therefore 

/,(a?)=8(4a?-2^), 
/ (a?) = 8» (2a?' - 6), 
/,(a?) = -8».20.a?, 
/,(a?)=-8M20. 

We see that the two last of the coefficients of the highest 
powers are negative. The roots therefore do not satisfy the con- 
ditions of stabiUty. 

As another example, take the equation 

/a? = a;« + aj* + 6a?* + 5a?' 4- lla?" + 6a? •{- 6. 

Here 

/j(a?) = a?*-6a?*+lla?'-6, 

/a(a?)=a?^-5a?»+6a?, 

/3(a?) = a?*-6a?^+6, 

/,(a?) = 0. 

Eeplacing /^ (a?) by the diflferential coefficient of/, (a?), we have 

/ (a?) = 4a?' -10a?, 
/,(a?)=10a?^-24a?, 

/(a?) =4.24. 

Here all the coefficients of the highest powers are positive, 
hence the roots satisfy the conditions of stability. 

It is clear that when the coefficients are numerical the rule 
given in Art. 5 is the most convenient, but when the coefficients 
are letters, the rule in Art. 9 will be found preferable. 

The process would be simplified by omitting the alternate posi- 
tive and negative signs of the terms in each line. 
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18. It may be interesting to express the two subsidiary func- 
tions fj (x) and f, (x) in terms of the roots of the given eqiuitian. 

Let ttj, a,, a^. . .a^ be the roots of the given equation/(aj) = 0, so 
that 

f{x) = {x-a^) (a;-aj(aj-a,)... 

Then it is evident that 
±f (a? V^) = a;* +p,a?*"* +p^otf"^ + . . . 

= J {os + a^)Qc + a^) ... + J(^-"«i) (^-^j) ••• 
It may be shown* that 

±A («^^^) ^p,x''-' + p^'^+ ... 

^,^ (a, + a)(a, + a)...(a, + aj ^ ^^^^^^ 

19. The foUowing propositions are not necessary to the main 
argument, but as they iUustrate geometrically the propositions in 
this chapter it has been considered proper to state them very 
briefly. The demonstrations wiU therefore be much curtailed. 

The equation being 

we put as in Art. 1, z^x-hy V— 1, and thus obtain two curves, 
whose equations expressed in polar co-ordinates are 

P -p^r^ cos w^ +;}jr"~^ cos (n - 1) ^ + . . . = 0] 
Q ^PJf'' sin«^ -k-p/"'^ sin (n - 1) ^ + . . . = Oj • 

These intersect in the radical points of the equation/ (^) =0. 

20. If we trace these curves we find that the curve P = 
has n asymptotes whose directions are given by cos nd = 0, i,e. 

^^ 1 TT 3 TT 6 TT 

n2' n2' n2' 

* Let us assume 

where A^, A^, &o, are constants whose values have to be found. Putting assa^, 
we have Pi^h*"^ +Pj^*"' + *o. = ii ^ (o^ - a J . . . (oi - a^). 

But since a^+pjpc^'^ + ... = (as - Oi) («- a,)...(a8 - aj, 

we have by putting x^ai and a^s - a^ 

Oi* -l>i«i**~^ + . . . = Saj (Oj + a,,) . . . (aj + a^). 

Subtracting the second of these results from the first we find il^ to have the value 
given in the text. . 
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These asymptotes all pass through the same point on the axis 

of ic. viz. a; = — — . It is also clear that only one branch of the 

curve can go to each end of an asymptote. Similar remarks apply 
to the curve ^ = 0, the directions of its asymptotes being given by 
sinw^ = 0, i.e. 

' n2' w2' n2' 

21. From these simple propositions we might, if it were 
worth while, deduce that every equation must have a root. The 
asymptotes of the two curves P = 0, Q = are alternatey and no 
two branches of the same curve can approach the same end of an 
asymptote. By sketching a figure, it may be easily shown that 
some branch of the F curve must cut some branch of the Q curve. 

22. Let us next consider the intersections of the curves P = 0, 
Q = 0. 

If we transform the origin to A, k, we put x = h + ^, y^k + rj. 

This is the same as expanding /(A + A? V— 1 +^ + ^V—l), and 
collecting into two parcels the real and imaginary terms. Let the 
expansion be 

where A^^A^, . . >are of the form 

c (cos a + sin a V— 1). 



If we put f + ^ V— 1 = r (cos 5 + sin ^ V— 1), we have 

P = CqCos o,^ + Cjf' cos {6 H- aj + c^ cos {20 + aj + ...' 
Q = Co sin a„ + c^r sin {0 + a^) + c^ sin {20 + a^ + ...} 

If the point (h, k) be a point of intersection c^ = 0. If the inter- 
section be a double point on either curve, the terms of the first 
degree must be zero, therefore c^ = 0, and the origin is therefore 
a double point on the other curve also. 

It is not difficult to show that if the intersection be a multiple 
point of any degree of multiplicity on one curve, it is a point of 
the same degree of multiplicity on the other curve. The tangents 
to these branches all make equal angles with each other, the 
tangents to the P and Q curves being alternate as we travel round 
the point of intersection. If the intersection be not a multiple 
point on either curve, the branches cut at right angles. 

/, Let us travel round a point of intersection along the circum- 
ference of a small circle whose centre is the point of intersection 
pi the direction in which is measured. Then it may be shown 
that as we pass from k P curve to a Q curve, P and Q have 
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opposite fnga^, and as we pass from a Q curve to aP curve the same 
sign. This is in fact merely Cauchy's rule for the changes of sign 

P 
of -jk when we travel round a radical point. 

23. Let us express the condition that there is no radical point 
on the positive side of the axis of y. This is the geometrical 
illustration of Art. 2. 

Draw a circle of infinite radius, and let it cut the asymptotes 
of the P curve in R, P,, P,,..,Pj, and the asymptotes of the 
Q curve in Q^, Q,,... Vj,. These pomta alternate with each other. 
Taking only those points which lie on the positive side of the 
axis of y, the P and Q curves may be said to begin at these points 
and are to intersect each other only on the negative side of the 
axis of y. The branches of the two curves must therefore remain 
alternate with each other throughout the apace on the positive 
side of the axis of y. Their points of intersection with the axis 
of y must be also alternate, and hence if we put a: = 0, in the 
equations P=0, Q = 0, and regard them as equations to find y 
the roots of each must separate the roots of the other. 

Conversely, we may show that if the intei-sections of the two 
branches are alternate on the axis of y, they cannot have intersected 
on that side of the axis of y on which the common intersection of 
all the asymptotes is not. This is the result arrived at in Art. 2. 

24. The following diagrams 'exhibit the forms of the curve 
P=0, Q = for a biquadratic. The dotted lines represent the 
asymptotes. 




A biquaJi'iitic witli four imagimry 
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A biqnadratio with equal imagiiuiy roots. 
25, [If jB=»+yV — 1, we have 

DifiTereutiating with respect, firstly to x and secondTy to y, we 



It easay tolbiis that ^ -4^ and 4^ = -I? ; 80 that hoth 
dm ay ax ay 

the functions P and Q satisfy the equation 

rf/ ^ die' "• 
The equation y(s) = gave us two curves which we have called 
P = and C =" 0. In the same way the derived equation /"(a) = 
will give U3 two other curves, which we may represent by P* = 
and Q' = 0. These we may call the derived P and Q curves. 

If as in Art. 22 we transform the origin to the point (A, k) we 
have 

P" = c, COS a, + 2(v cos (^ + 0^ + 3c,r* cos (25 + a,) + . . . 
Q' = Ci sin Oi + 2c^ sin (fl + aj + 3c/ sm{20 + a^ + ... 
If the origin he at a point of intersection of the curves P= 0, 
P* = which is not a double point on the first of these curves, we 
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have Co cos a© = and cos a^ = 0. Hence the tangent to the curve 
P = at the point of intersection is parallel to the axis of x. 
Conversely if we move the origin to any point on the curve P = 
at which the tangent is parallel to the axis of x, we find that the 
curve P' =0 passes through the origin. Hence the derived P 
curve passes through all those points on the curve P = at which 
the tangent is parallel to the axis of x, and all those points on the 
curve ^ = at which the tangent is parallel to the axis of y. In 
the same way the derived Q curve passes through those points on 
the curve Q = at which the tangent is parallel to the axis of x, 
and those points on the curve P=0 at which the tangent is 
parallel to the axis of y. 

If CoCosao=0, and 0^=0, the origin is a double point on the 
curve P=0 and the origin also lies on both the curves P' = 0, 
Q' = 0. So that both the derived curves cut the curves P = and 
Q = in their double points. In other words, these double points 
are radical points for the derived equation /'(«) =0. In the same 
way, any multiple point on either curve is a multiple point of one 
degree less multipUdty on both the derived curves. 

If a finite straight line AB be drawn parallel to the axis of x 

joining two points A, B on the same or on diflferent branches of 

the curve P = 0; this finite straight line must cut one or more 

branches of the derived P curve. For it is clear that if P vanishes 

rfP 
at A and P, P which i& equal to -r- cannot keep one sign between 

A and B, and must therefore vanish somewhere between A and P. 
If A and B be adjacent points, Le, if there be no other points 
between A and B belonging to the curve P = 0, then the straight 
line AB must cut an odd number of branches of the derived P 
curve. In the same way if a straight line CD be drawn parallel 
to the axis of y joining two adjacent points on a derived Q curve, 
this straight line must also cut an odd number of branches of a 
derived P curve between G and D. 

By considering a tangent as the limit of secant, it again follows 
that if a tangent be drawn to the curve P=0 parallel to the axis 
of X, the derived P curve must pass through the point of contact. 

Let P be a radical point on the derived curve /'(«) = 0, and let 
it not be a double point on either of the curves P= 0, Q = 0. Let 
a straight line be drawn from R in any direction cutting the 
branches of either of the curves P=: 0, (? = in the points A^, A^ 
&c. Then we may show that 

so that the polar line of R with regard to either of the curves 
P = Q, Q = is at infinity. 
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The positions of the radical points of the derived equation 

jT («) = relatively to any branch or branches of the curves P = 0, 

Q = may be found by the use of Cauchy's theorem. If the 

point spoken of in Art. 1 travel along a branch of the curve P=0, 

it is easy to see that 77 = ;; • If i^ travel along a branch of the 

F dx 
curve ^ = we have 77 = ~ ^ • If ^&tcl any contour be partly 

bounded by branches of these curves, the simplest inspection of 
the points at which the tangents are parallel to the axis will 

F 

determine the changes of sign of yy as it passes through zero. 

If another part of the contour be an arc of a circle of infinite radius 
whose centre is the origin, the changes of sign through zero will 
be from + to — and their number will be indicated by the number 
of asymptotes of the derived P curve which cut the arc] 

26. The use of Watt's Governor in the steam engine is too 
well known to need description. It has however, as commonly 
used, a great defect. It is sometimes of importance that the 
engine should continue to work at the same rate notwithstanding 
great changes in the* resistances. Suppose the load suddenly 
diminished, the engine works quicker, the balls diverging cut off 
the steam, and the engine, after a time, again works uniformly, 
but at a different rate from hefbre. The balls as they open out 
or close in are usually made to describe circles. Let them now be 
constrained to describe some other curve which we may afterwards 
choose so as to correct the above defect. If this curve be a parabola 
and the balls be treated as particles, it is clear from very elemen- 
tary considerations that these will be in relative equilibrium only 
when the engine works at a given rate. This principle is due 
to Huyghens, see Astronomical Notices, December, 1875. It is 
now proposed to determine the condition of stable oscillation about 
a state of steady motion. 

Two equal rods AB^ AE are attached at A by hinges to a 
small ring which can slide smoothly along a vertical axis. The 
ring is attached by a rod to the valve and can thus govern the 
amount of steam admitted. Two equal balls are attached at B 
and B, and the centre of gravity G of the rod AB and the ball B 
is constrained to describe some curve. To represent the inertia of 
the engine we shall suppose a horizontal fly-wheel attached to the 
vertical axis whose moment of inertia about the axis is /. Let the 
excess of the action of the steam over the resistance of the load 
be represented by some couple whose moment about the vertical 
axis is f{6)y where 6 is the inclination of the rod AG to the 
vertical, and / is a function which depends on the construction of 
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the engine. Since the steam is cut off when the balls open out, 
it is clear that '^^ ^ is negative. 

There may be also some resistances which vary with the • 

velocity. Let these be represented by a couple B -j tending to 

retard the motion round the vertical axis and a couple round A in 

JO 

the plane BAB equal to mC -r- ^ where ^ is the angle the 

vertical plane BAB makes with a fixed vertical plane. 

Let m be the mass of either sphere and rod; k the radius of 
gyration about an axis through perpendicular to the rod, and k' 
that about the rod, let i = -4 O, Then the equation of angular 
momentum ^ves 

Let the steady motion be given by ^ « a, -~ =^n, and let the 

oscillations be represented by ^ = a + a?, ^^n-\-y. We have then 

m=f{a)+f{a)x, 
/(a)=5n. 
The equation then reduces to 

{/ + 2m (FT? sin' a + T cos' a)} ^ + 2m (f + A;' - jfc") sin ian ~ 

^f{a)x-By. 
This equation may be briefly written 

Let z be the altitude of Q above some fixed horizontal plane. 
Then if T be the semi vis viva 

2 r= /^'^ + 2m {(A:' + Z") sin»^ + A?" cos»5} f^'^ + 2mFd'" 

+ 2m|(^^y+Pcos'^U'«, 

If U be the force function, omitting the couple of resistance, 
we have U— — ^mgz. The virtual moment of the couples of 
resisfence being 2mC0'Bd, the Lagrangian equation of motion 

becomes 

d dT dT_dU ^^ 

dtWdd^de"'^'^^^' 
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Substituting for Tand J7we have when the motion is steady 



9 

^..^(P + A;»-r)sin*^ = ;5. 

Since ?sin ^ is the distance of Q from the vertical axis, we see 
that the path of O must be a parabola. The semi latus rectum is 

which, we notice, is independent of the radius of the balls. The 
length of this latus rectum must of course be adjusted to suit the 
particular rate at which the engine is intended to work. 

When the system is oscillating about the state of steady motion 
we have, putting 

a» = P + ^(Z» + A;'-Jfc'7sin*a, 
9\ 

and rejecting the squares of x and y, 

(ft' + a'cos'a) J + cJ-(P + ft*-A;'Osin2any = 0. 

The term a?, it will be noticed, has disappeared from the 
equation. This equation may be brieiSy written in the form 

Eliminating y from the equations of motion we have 
AH^^-{AC+BH)^ + {BG^-EL)^^-^FLx = 0. 

The coefficients are all positive, the necessary and sufficient 
condition of stability is therefore 

{AC^ BH){BG-^EL) > AFHL. 

In some clocks to which Watt's Governor is applied, there is a 
special arrangement which causes (7 to be much greater than B. 
See the Astronomical Ifotea, xi„ 1851, and the Memoirs of the 
Astronomical Society, Vol. xx. Neglecting therefore B, we have 

CE>FH. 

:. 2(7mn(P + jfc'-ifc'^sin22>J^(i!^ + a»cos»a). 
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Fommtion of the eqimtiona of steady motion and of smoM oscillation 
where Lagrang^s method may he used. Arts, 1 — 5. 

The equations being all linear the conditions of stability are caressed 
by the character of the roots of a determina^tal eqiuition of an even order. 
Art, 6. 

Mode of easpa/nding the determinant. Art, 7, 

A method of Jmding the proper coordinates to make the coefficients of 
the Lagrangia/n function consta/nt. Arts, 8 — 10. 

How the Harmonic oscillations ahout steady motion differ from those 
about a position of equilibrium. The forces which cause the difference are 
of the nature of centrifugal forces produced by an im^aginary rotation 
about a fixed straight line. Arts, 11 — 19, 

Reduction of the funda/mental determinant to one of fewer rows by the 
elimination of aU co-ordinates which do not appear except as differential 
coefficients in the Lagrangian function; vrith an example. Arts, 20 — 23. 

Formation of the eqtuitions of Motion and of the determinant when 
the geometrical eqtuitions contain differential coefficients, so that Lagrange's 
me^iod comnot be used; with an example. Arts, 24 — 27. 

1. Let the system be referred to any co-ordinates f , 17, f, &c. 
The general expression for the kinetic energy is 

where P, Q, &c. are known functions of f, 17, f, &c. and accents 
have their usual meaning. Let us suppose the system to have 
some motion represented by 

and when disturbed, we wish to find the oscillations about this 
motion. To effect this, we put 
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where ^, ^, &c. are all small quantities. Substituting and expand- 
ing T in powers of ^, 0, &c., we find ' 

+ ^{A^0' + 2AJ<I>+...) 

+ &C. 

In the same way we may make an expansion for the Potential 
Energy of the forces, viz. 

^^{E,^0'^2EJj> + .„), 

when these two functions are given the whole dynamical system 
and the forces are known; and we may form the equations of 
motion by Lagrange's method. 

2. We shall here however limit the question by supposing 
that the motion about which the system is oscillating is what has 
been called in Chap. I. steady. The analytical peculiarity of such 
a motion is that when referred to proper co-ordinates, every coeflS- 
cient in each of these two series is constant, i.e. independent of t 
As already explained the physical peculiarities are that the vis 
viva is constant throughout the steady motion and the same oscil- 
lations follow from the same disturbance at whatever instant it 
may be applied to the motion. A method of discovering the 
proper co-ordinates, if unknown, will be given a little further on. 

3. In order to form the equations of motion we must now 
substitute in Lagrange's equations 

dtdO' d0^d0^ ' 

&c. = 0, 

rejecting all the squares of small quantities. The steady motion 
being given by 0, 0, &c. all zero, each of these must be satisfied 
when we omit the terms containing 0, (f), &c. We thus obtain the 
equations of steady motion, viz. 

A, = E„ 

&c. = &c. 
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These equations may be simply formed in any case by the 
following rule. Putting i = T— F, so that L is the difiFerence 
between the kinetic and potential energies, expand L in powers of 
the co-ordinates 0y (f>, &c. regarding ff, ^\ &c. as zero. The required 
relations are obtained by equating the coefficients of the first 
powers to zero. This rule may be also expressed thus. Let L he 
the general expression for the excess of the kinetic energy over the 
potential energy of a dynamical system in terms of its n co-ordinates 
f, 17, &c. Let this system be moving in steady motion with constant 

valms ?/" jf > 7^ > &c.* Then substituting these constant values in 

the general expression for L, the relatione between the constants of 
steady motion are given by 

dL ^ dL ^ « . 
jTj = 0, 3- = 0, &c. 
df d77 

In this way we obtain in general as many equations as there 
are co-ordinates. Usually the coefficients in the expression for T 
are constant because some of the co-ordinates are constant in the 
state of steady motion, and the other co-ordinates appear in the 
expressions for T and V only as dififerential coefficients. In such 
cases we have clearly fewer equations than co-ordinates to con- 
nect the constants of steady motion. We have then a system of 
possible steady motions which we may conveniently term parallel 
steady motions^ 

4. To obtain the equations to the oscillatory motion, we retain 
the first powers of 0, ^, &a We thus obtain a series of equations 
of which the following is a specimen : 



+ 



5.S+(c;.-cj|-^,+^Jt 



df • ^^«^ ^^^ dt 
+ =0. 

To solve these we write 

e^M^e"^, <f> = M,e'^, '^^M^e'^\ &c. 

Substituting we obtain on eliminating the ratios M^\ M^\ M^ &c. 
a determinantal equation, viz. 

* Since we may in Art. 1 change the co-ordinates from f, 17, f, &c. to f j, ly^, fj, 
Ac, where ^=/(^j) 7; = 1^(7; 1) &c., it is clear that the steady motion can be always 
expressed by constant values of the differential coefficients of the co-ordinates. 
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By^m^-A^ + E^ 


B^m'-A^ + E^ 


Bj^m* -Aa + Eji 


&c. 


B^m*-A^ + E^ 


B^m'-A^+E^ 


B„m*-A„+E„ 
+ (C„-(?„)m 


&c. 


B^m'-A,,+E„ 


B„m!'-A^+E, 


B„m'-A^+E^ 


&c. 


&c. 


&C. 


&c. 


&c. 



= 0. 



This equation will be referred to as the Determinantal equa- 
tion. 

6. If we refer to the -equation formed by this determinant 
and read it in horizontal lines, we have of course the several 
equations of motion, each term being the coefficient of 0, ^, >|r, &c. 
in order. In this form the equations may be reproduced by the 
following easy rule. 

Taking the expression for T— V as given in Art. 1, let us 
consider only the terms of the second order, those of the first 
having been already used to determine the steady motion as ex- 
plained in Art. 3. Separate from the rest, the even powers of 
e, 0, (f>, (f>\ &c. and write for ^^ ^^', &c. - 1^6^ - i)'^^^, &c., so 

that D will stand either for -rr or for the m in the determinant 

at 

of Art. 4, when we write = M^e"^, <j> = M^e*^, &c. Let the sum 

of these terms be called P, so that 

P = J M„ - E^, - B,,I>') 0' + {A,, - E,, - B,,If) 64, + &c. 

Let the remaining portion of the terms of jT— F, viz. those con- 
taining both 6, (f>, &c. and 0", <!>', &c.,-be called Q, so that 

Then the several equations may be formed from the rule 

dtf^ d0 ^ dff ~ 

d<f> d(f> dcf} 

&c. = - 

In applying this rule no accented letters will occur except in 
the second term of each equation. If we wish D to stand for m 



r . 
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in the determinant, we must regard ff, f . &c. as abbreviations for 
Ddy D(f>, &c. If we wish to use the equations themselves, we 
replace 1)0, Dtf), &c. by 0', <j>, &c. 

[The determinant may also be found by another rule. Taking 
as before only the terms of the second order in the Lagrangian 
function X = T— F, let us separate the terms of the form 

Q = CJff + CJ<I>' + C^ff<f> + &c. 

In the remaining terms put ^ = ^m V— 1, <^' = <f>m V— 1, and 
so on, and write down the discriminant If the system oscillates 
about a position of equilibrium the terms represented by Q are 
absent and the discriminant thus formed will be the determinantal 
equation giving the required values of m. But if the system 
oscillate about a state of steady motion we must modify the dis- 
criminant by adding some quantity derived from Q to each term. 
To find this, write above the columns ff, (f), &c. and before the 
rows 0, (jj, &c. Consider any term, say the term in the column 
0' and row 0. We must add to that term ( (7,^ — (7,^) m, where 
U,j — Cj^ is the excess of the coefficient of (f>'0 above the coefficient 
of <f>0' m the expression for Q. Since the determinant is un- 
changed by writing — m for m, we may, if preferred, add the 
excess of the coefficient of ff(l> above the coefficient of 0(f> provided 
we adhere to one order throughout] 

6. If in the determinantal equation we write — m for m, the 
rows of the new determinant are the same as the columns of the 
old, so that the determinant is unaltered. When expanded, we 
shall have an equation which contains only even powers of m. 

The condition of dynamical stability is that the roots of this 
equation should be all of the form 

If we write m*=— p^ the roots of the transformed equation 
must be all* real 

• In the case of equal roots of the. form m= + ^V— 1, it has 
been shown in Art. 5 of Chap. I. that it is necessary for stability 
that the proper number of minors in this determinant should 
vanish. If there be two equal roots, these roots must make all 
the first minors zero ; if three equal roots, all the first and second 
minors must vanish, and so on. 

7. When the system depends on many co-ordinates the labour 
of expanding this determinant is often considerable. Methods of 
evading this in certain cases will be given in the following chap- 

R. A. 4 
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ters. But when the development is necessary we may proceed in 
the following manner. Let the determinant be written 



B^^rri^ + a^^ 



&c. 



and let 



+ ^12^ 

B^rn^ + «„ 
&c. 



^ 



+ ^98^ 



&c. 



-5= -Bjj 2 + B^0<f> 4- ..., 



We know that the detenninant when expanded is of an even 
order, hence all odd powers of m must finally vanish. Let us 
expand the determinant in powers of the F'b, The first term is 
the discriminant of Bw? + a, this term is independent of the ^s. 
The terms which contain the first powers of the J^*s are obtained 
by erasing any one line of this discriminant and replacing it by 
the corresponding F terms. But these terms all vanish and we 
need not describe them minutely. The terms containing the 
products and squares of the -F's may be obtained by erasing every 
two lines of the discriminant and replacing them by the corre- 
sponding i'^ terms. Thus if we erase the two first lines we have 
the determinant 






F,,m 


■^u^W 


&o. 


-F„m 





F^m 


&c. 


B^m' + a„ 


5««»'+o„ 


B„m* + a„ 


&c. 


&c. 


&c. 


&c. 





and so on for all the other rows taken two and two. The terms 
which contain the cubes and all odd powers of the JPs vanish, 
while the terms which contain the fourth powers may be obtained 
by. erasing four lines of the discriminant and replacing them by 
the corresponding JPs. 

When the determinant has been expanded, we have an equa- 
tion of an even order to find the values of m. We may therefore 
employ the short method of Art. 5, Chap, in., to obtain the 
Sturmian functions. 

8. Necessary and sufficient tests of the stability of the motion 
of a system of bodies are given in the preceding pages. But it 
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is assumed, as explained in Art. 2, that the co-ordinates have been 
properly chosen. They are supposed to have been so chosen that 
the coefficients in the expanded Lagrangian function are all 
constants. When this is not the case we must discover the 
proper co-ordinates to which the system must be referred before 
we can apply the test of stability. But when the motion is steady 
this is not difficult. 

There are obviously many such systems of co-ordinates, and 
one set may generally be found by a simple examination of the 
steady motion. If there are any quantities which are constant 
during the steady motion> they will often serve for some of the 
co-ordinates. Others may be found by considering what quan- 
tities appear only as differential coefficients or velocities. Practi- 
cally these will be the most convenient methods of discovering 
proper co-ordinates, since no further change will then be necessary 
and we may at once form the determinant of stability. But if 
these methods fail we may adopt the following analytical method 
of transforming (where possible) the general Lagrangian function 
with variable coefficients into one with constant coefficients. 

* 

9. Let the Lagrangian function be 

i = i, + ^,^ + A<^ + &c. + Bjfl' + -B,f + &c. 

+ J5,/' + 5,/f + ... 

where the coefficients are all functions of t and the co-ordinates 
By 0, &c. have been so chosen as to vanish along the steady motion. 
We have therefore for the steady motion 

&c. =0. 

The oscillations about the steady motion are given by the 
terms of the second order. Our present object is to transform 
these to others with constant coefficients by the following substi- 
tutions : 

-p^x +p^ ^p^z + &c., 

(f^^q^x + q^ + q^ + Scc, 

&c. = &c., 

where the p's, j's, &c. are functions of the time at our disposal 

Substituting and equating the coefficients of x'\ y", &c. to 
unity, we have as many equations of the form 

iB,^' + iiB„^ + B„pq + ... = l... (1) 

4—2 
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as there are co-ordinates. Equating the coefficients of the pro-» 
ducts a?y , x'z\ &c. to zero we get n — 5— equations of the form 

Ail^i-P« + -S«M.+-Si»(M2+i'a?i) + -M =0 (2), 

supposing that there are n co-ordinates. 

Equating to the constants a^, a,, ... the coefficients oixx\ yy\ 
&c, having subtracted the differential coefficients of (1) we have 
n equations of the form 

(Cn-iA.')i''+(C„-4-B„')2'+(^«+C«-A.')i'2+-=«-(3). 
Adding the coefficients of xy and xy and subtracting the 
differential coefficients of (2) we have n —^ — equations of the 
form 

Equations (1), (2) and (4) give n* equations to find the n* 
quantities p^p^^ &c., q^q^ &c. The solution of these equations is a 
purely geometrical prolblem. If we construct the two quadrics 

and refer them to their common conjugate diameters, by writing 

0^PiX+pj/+... 

(l> = qiX + qj/-\' ... 

&c. = &c., 

making the first quadric to become what we may call a sphere by 
projection; the values otp^p^y &c., q^q^ &c. thus found are the values 
required to make some of the coefficients in the Lagrangian 
function become constant. These values must of course make all 
the other coefficients of the second order in the Lagrangian function 
constants also, and thus we havew(w+]) analytical conditions 
that the motion should be steady. 

It might be supposed that greater generality would be obtained 
by replacing the zero's of equations (2) and (4) or the unities of 
(1) by arbitrary constants. This may be oonvenient in practice, 
but as we know that by a subsequent real change with constant 
values of jpip^, &c., q^q^, &c., we can render them zero or unity, it 
simplifies the argument to perform the two transformations at 
once. 
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10. The geometrical problem just alluded to admits of a real 
solution whenever one quadric can be projected by a real projection 
into a sphere. The problem then becomes that of finding the 
principal axes of the other. This is just our case, since the ex- 
pression 

is necessarily positive for all values of ^^'. 

It is unnecessary to describe here the mode of solving this 
problem. It is sufficient to say that it may be reduced to the 
solution of the symmetrical determinantal equation 

■^11 "* f^-L^ny -^M — ^-^12 • • • 



•••••• 



= 0, 



where D^, D^, &c. are the coefficients of 0*, 6(f>^ &c. in the second 
of the quadrics. The roots of this equation are known to be real 
when the suppositions just mentioned are satisfied, 

il. In order to examine the fundamental determinant in 
Art. 4 a little more closely, let us suppose it reduced to depend on 
three co-ordinates. We may then have the advantage of a geo- 
metrical analogy. Let the co-ordinates be f 17, f, and let the 
equations of motion be written 







<? . . x. <p 



( .4 r\ .4 j'^(^-^-^")'-''- 

Let a geometrical point P move in space so that its co-or- 
dinates referred to any axes are f, 17, f. Then the position and 
motion of the point exactly give us the position and motion of 
the system. 

12, Looking at the equations of motion just written down we 
see that they are similar to those which give the oscillations of a 
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system about a position of equilibrium, but that there are in 

addition terms E -^, F-^, &c. The general effect of these termSi 

as wiU appear from what foUows in the subsequent chapter, is to 
increase the stability. If we transpose these terms to the other 
sides of the equations, we may regard them as impressed forces 
acting on the system, whose resolved parts in the directions of the 
axes f , 17, ^, are 

at at 



at at J 



We see at oace that 



EX+FY+GZ=0 





at at at 



so that these forces are at once orthogonal to the path of th^ 
representative point P and also orthogonal to the straight line 
whose direction cosines are proportional to E, F, 0, These forces 
are therefore of the nature of centrifugal forces^ as if they were 
produced by the rotation of the system about this straight line. 

13. We may show that the straight line {E, F, Q) is fixed in 
space. To prove this, let us transform our co-ordinates from ^,17, f 
to a?, y, Zy where a?, y, z are connected with f , 17, { by any linear 
relations, such as 

^==a^ + hj/ + c^, 

Let the portion of the Lagrangian function under consideration 
(Art. 1) be 

then (? = (7,, - C;„ E^G^^C^,F^ G,, - C,,. 

Substituting for f , 1;, f their values in terms of a;, y, z, we find that 
the difference between the coeflScients of xy' and x'y is 



«»«. 


+ E 


«««8 


^■F 


«»Ol 


Kh 




h\ 




hh 
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with Bimilar equations for E and W. But if /t be the determi- 
nant of transformation 

/*« = ? 



a,o, 


+ ? 


«.«8 


+1? 


a.«i 


hK 




\K 




J, &. 



with similar equations for x and y. The ratios of -E?, jP, Q are 
therefore transformed as if they were co-ordinates. If the trans- 
formation be a real transformation of Cartesian co-ordinates, let 
lengths each equal to unity be measured from the origin along the 
axes (?f, 0% 0£ thus forming a tetrahedron whose volume is F. 
Let a similar construction be made for the new axes, forming a 

V 
tetrahedron of volume Y, Then* fi^-r=z. Hence the quantities 

E F G- 

T^, T^, ^ may be transformed as if they were lengths measured 

W W C" 
along the axes and become y7 , y? , y? . If both systems of co- 
ordinates are rectangular we have F= ^ , F = ^ . 

w jr Q. 
Let ft) he the resultant of y^^ y, p., then co may be regarded 

as a fijxed length measured from the origin along a straight line 
fixed in space. 

Let V be the velocity of the representative particle, the angle 
between the direction of this velocity and the axis whose direction 
cosines are proportional to E, F, O, Then the resultant of the 
forces X, F, Z \& easily seen to be 2v V(o sin acting perpendicular 
to the axis and to the dii'ection of the motion. We might call the 
straight line {EFG) the axis of the centrifugal forces, 

[* Let (^1 T/i fi), ftj 17, fa), (i. iji fg) be the co-ordinates of three points A, -B, C re- 
ferred to any oblique co-ordinates. Let us find the yolume V of the tetrahedron 
of which these and the origin are the angular points. Since the volome vanishes 
when any angular point as C Hes in the plane containing the origin and the other 
two A, By the expression for the volume must contain the factor 



/*= 



^1 ^s ^8 

I?! l/t "tz 

ti ta ^8 



The volume is evidently an integral rational function of the co-ordinates when 
the axes are rectangular and the plane AOB is taken as the plane of xy, it 
easily follows that this is true for aU axes. Since this function cannot be of more 
than the third order, we have V'=Mfiy where M is independent of the co-ordinates 
oiA,B,C, When the points A^ B, G are on the axes at unit distances from the 
origin, let V be the volume of the tetrahedron. Li this case Ab=l, and ,\ M^V, 
We have therefore in all cases V = 7/a. 

Li the text, let the extremities of the unit lengths measured along the axes of 
90, y, z be called A^ B, C. Then the (^vi) co-ordinates of A, By C are (a,, a,, a,), 
(61, 6,, 6,), (Ci, c,, c,), respectively. Hence by what has just been said v'=Vfi.] 
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14. The expressions for the co-ordinates in terms of the time 
will in general contain as many periodic functions as there are 
co-ordinates. If the initial conditions are such that each contains 
one and the same periodic function, the motion recurs after a con- 
stant interval and the system is said ta be performing a simple or 
harmonic oscillation. 

If the system be oscillating about a position of equilibrium, 
with a Lagrangian function 

we know* that the harmonic oscillations are represented by rectili- 
near motions of the representative particle, and that these are 
along the conamon conjugate diameters of the two quadrics ^^^ 
equations are 

-4u| + ^i2f»7 + •••=« 

where a and 5 are two constants chosen to make the quadrics real. 
Let us consider what are the harmonic paths of the representative 
point when the system is oscillating about a state of steady motion. 

In any harmonic vibration we have f = i.cos (X^ + a) with 
similar equations for 97 and f. Hence 

df'"^^' w — ^^' w — ^^- 

Substitute these in the equations of Art. 11. Dififerentiate and 
substitute again. Multiply by |, 77, f, add and integrate, we obtain 

where c is some constant. The harmonic path lies on this quadric, 
which has a common set of conjugate diameters with the two 
quadrics a and 6. 

If we resume the result of the substitution of -^ , &c. in the 

air 

equations of Art. 11, and multiply by j&, -F, respectively and add, 

we obtain 

[{B^E^ B^+B^G)^ + &c.] V + [{A^E+ ...) f + &c.] = 0, 

* A short paragraph in Thomson and Tait*s Natural Philosophy y page 273, is 
thinly notice of this which the author has discovered. 
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which is a plane, and is diametral to the straight line [EFO) with 
regard to the quadric c. 

The harmonic paths are therefore ellipses. The three harmonic 
planes are diametral to the same straight line and this straight 
line is fixed in space, being the axis of the centrifugal forces. 

If we eliminate \ between the equations to the plane and the 
quadric c, we get a cubic surface on which the three harmonic 
conies lie. 

If Ey Fy are zero, which is the case when the system oscillates 
about a position of equilibrium, the quadric c- becomes a cylinder. 
This may be conveniently shown by referring the system to such 
co-ordinates that the coefficients B^y B^y B^y Aj^y A^, A^ are all 
zero. In this case the diametral plane of every straight line passes 
through the axis of the cylinder. The harmonic oscillations are 
therefore rectilinear. 

If B be the length of that semidiameter of the quadric (c) 
which is parallel to the fixed straight line {E, F, Q), it may be 
shown that the 

Product of the axes ] R 2c 

of the quadric c ) JE^ + i^" + G* ' ^ * 

If Ey Fy a are all zero, and their ratio is indeterminate, B is 
any diameter. Hence one of the axes of the quadric (c) must be 
infinite and the quadric will be a cylinder. 

[If the quadric (c) be a cylinder and Ey F, G are not all zero, 
we must have either \ zero or B infinite. In the latter case the 
axis of the cylinder will coincide with the axis of the centrifugal 
forces.] 

The quadric (c) has also the following geometrical property. 
Let the lengths of semidiameters of the quadrics (a) and (5) drawn 
parallel to the axis of the centrifugal forces be p and p. Through 
the intersection of these quadrics describe a quadric so that the 



Product of 1 2Va5 Vp"* pV 

^/WTT'VCP' Al _lVi. -1\ ' 

Vvp'" ipAi? pV 



}= 



its axes 

,_ . _ . ^ A/ I 

This quadric is similar to the quadric (c). 



13. The introduction of the representative point to exhibit 
the motion of a system may appear somewhat artificial. If how- 
ever we properly choose the co-ordinates the particle moves exactly 
as a free particle, and we might reduce the problem of finding the 
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oscillations of a system to a problem in Dynamics of a particle. 
Refer the quadric 

to its principal axes and let the equation thus changed be 

Since -B^ , B^, &c. are positive quantities, we may put 

The quadric has thus been "projected** into a sphere. Let a?, y, z 
be now chosen as the co-ordinates of the system and let the 
Lagrangian function be expressed in the form 

the terms of the first degree being omitted as not necessary to our 
present purpose. The three equations of motion at the banning 
of Art. 11 take the form 

&c. = 0, 

which are the three equations of motion of a free single particle 
whose co-ordinates are x, y, z under the action of forces whose force 
function is 

and a force acting perpendicular to the path and also perpendicu- 
lar to a fixed straight line, the force oeing proportional to the 
velocity. 

16. As an illustration of this theory, let us here make a short 
digression. However the particles of light may oscillate, whether 
in a rotatory or linear manner, we know the motion is related to a 
certain plane called the plane of polarization. It may be shown 
that any harmonic oscillation about a position of equilibrium may 
be represented by a rectilinear oscillation of the representative 
particle. Let us represent the motion at any point of the ether 
by a rectilinear oscillation in a direction perpendicular to the 
plane of polarization. This would be Fresnel's Vibration. The 
representative particle, as just shown, wDuld not necessarily move 
as if it were a free single particle. But let us assume (and a proof 
is not necessary to our present purpose) that when the oscillation 
is drawn as above described the motion in the plane of the front 
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is the same as that of a free particle, while that perpendicular is 
not free. On this assumption we see that Fresnel in his theory 
of double refraction is justified in taking actual instead of relative 
displacements, for it is the representative particle he is considering. 
He also neglects the force normal to the front, for the particle 
moves as a free particle only in the plane of the front. These 
general remarks are not meant to explain FresneVs theory, but 
merely to show how the representative particle may be used to 
replace a complicated motion. 

17. [When a system is performing a harmonic oscillation 
about a state of steady motion or about a position of equilibrium, 
the motion repeats itself continually at a constant period, that is 
to say, the values of the co-ordinates recur at this interval. This 
is the chief peculiarity of a harmonic oscillation. When the 
oscillation is about a position of equilibrium, the representative 
particle oscillates in a straight line whose middle point repre- 
sents the position of equilibrium. Thus the system passes through 
the position of equilibrium twice in each complete oscillation. 
When the oscillation is about a state of steady motion the path 
of the representative particle is an ellipse whose centre is at the 
point occupied by the system in steady motion at the same in- 
stant. Thus the system does not in general ever coincide with 
the simultaneous position of the system in the undisturbed or 
steady motion. When a system is disturbed by a small impulse 
from a state of steady motion, it will in general describe a com- 
pound oscillation made up of at least two harmonic oscillations, 
at the instant of disturbance these two neutralize each other so 
that in the disturbed and undisturbed motions two simultaneous 
positions are coincident. But it is clear that this cannot occur 
again unless either the periods of the two harmonics are com- 
mensurable or the period of one of them is infinite.] 

18. [In some cases the ellipse degenerates into a straight line. 
Thus if the quadric (c) be a cylinder the diametral plane of the 
axis of the centrifugal forces will pass through the axis of the 
cylinder, and thus the harmonic oscillation corresponding to this 
particular value of \ will be rectilinear. In this case the system 
twice in each oscillation passes through the position it would have 
occupied at the same instant in the undisturbed motion. 

The quadric (c) has a common set of conjugate diameters with 
the quadrics (a) and (i). Hence if (c) be a cylinder, its axis must 
be parallel to one of the three common conjugate diameters of (a) 
and (i). If we refer the quadrics (a) and (i) to their common 
conjugate diameters, they take the form 
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The cylinder which passes through their intersection and has 
its axis parallel to the diameter fis found by eliminating f* between 
these equations. We see therefore that -BggV + A J = 0. If then 
the axis of the cylinder cut the quadrics (a) and (6) in D and jy 
respectively, we find that for this oscillation 

,_ od;^ -a 

" OD^ ' b ' 

It has already been shown that wh^n this value is finite, the 
direction of ODD' is along the axis of the centrifugal forces.] 

19. [In some cases two or more of the values of X are zero. 
In these cases the co-ordinates will have terms of the form nt + e, 
where n and e are two small constants. When, as explained in 
Art. 3 of this Chapter, there are several parallel states of steady 
motion, these terms imply that the motion is stable about a state 
of steady motion very nearly the same as the undisturbed motion 
but not coincident with it. The actual undisturbed motion, unless 
n is zero, is unstable in the sense that if a proper disturbance 
be given to the system, the system will depart widely from the 
positions it would have simultaneously occupied in the undisturbed 
motion.] 

20. In many cases of small oscillations it will be found that 
the Lagrangian function T^V is not a function of some of the 
co-ordinates as 0, <f>, &c. though it is a function of their dififerential 
coeflBcients ff^ </>', &c. In such cases the steady motion will be 
usually given by constant values of these dififerential coefficients, 
while the other co-ordinates as f, 77, &c. are also constant. It is 
evident that the determinantal equation of Art. 4 is needlessly 
complicated. It is clear that there will be as many pairs of roots 
equal to zero as there are co-ordinates 6, 0, &c. It will be an 
advantage to eliminate 6\ ^', &c. altogether from the Lagrangian 
function, and to find the remaining roots by operating only with 
the co-ordinates f , 17, &c. We shall thus obtain a determinant with 
just as many rows as there are co-ordinates of the kind f , 77, &c. 

Let ij be the Lagrangian function expressed as a function of 
ff<}>y &c. ^, ^V> &c. Let L^ be its value when ff^' are elimi- 
nated, so that iy, is a function of ^, ^rj', &c. only. To effect this 
elimination we have the integrals 

dT_ dT _ o 

where c^, c,, &c. are constants. Then 

dL,_dL dl^dl dL,d£ 
d^~d^^ dd' 'd^^d^' dt 
_ di, dff d^' 



N 
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dL^_dL dff d<i>' 



«?f d^ ^ ' df ' df 



d^dL^ dL,_d^dI^ dL /d^dff_d9\ 
dt df df ~ dt d^' df "*■ "• U< df rff/ 



But ^^_^t_0 

^"* d< df df - "• 

Hence if we take 

y = 'L-c.0 -cjh'-Jkc 

and eliminate ff, ^' 6y AeZp o/" ^c integrals j^ = Cp j^i/^^a* 

w;6 may ^rca^ L' just as we do the Lagrangian function L. The 
equations giving the smaU oscillations about ihe steady motion 
tvill be 

d dL' dL'__ 

dtdr"""dF""''^^-"^' 

The function L' may be called the modified Lagrangian function, 

[It should be noticed that this is equivalent to a partial use of 
Hamilton's transformation of Lagrange's equations. Sir W. R. 
Hamilton eliminates all the dififerential coeflScients ^, 0', &c, by 

the help of equations of the form jff = % ^a' ~^» ^^' ^^^^^ ^> v> 
&c. are made to be new variables*. In our transformation only 

* The Hamiltonian transfonnation of Lagrange's equations bears a remarkable 
analogy to the transformation of Beciprocation in Geometry. This may be shown 
in the following manner. 

When the system has three co-ordinates 0^ </>, ^, we may regard 0^, 4>\ r// as the 
Cartesian co-ordinates of a representative point P. The position and path of P wiU 
exhibit to the eye and will determine the motion of the system. Let u, v, u; be the 
Hamiltonian variables, so that 

dT. dT. dT. 

''=dP' "=50- "'=#" 

where Tj is the semi vis viva expressed as a function of 0, 0, yp, 0', </>', }//, Then 
u, V, w may be regarded as the co-ordinates of another point Q whose position and 
path will ^80 determine the motion of the system. 

If the semi vis viva be given by the general expression 

it is clear that the point P always lies on the quadric T^ = 17 where U is the force 
function and the co-ordinates 0, <t>, \f/ have their instantaneous values. The point 
Q must therefore lie on another quadric which is the polar reciprocal of the first 
with regard to a sphere whose centre is at the origin and whose radius is equal to 
»J2U. The equation to the reciprocal quadric is therefore 

tt V w =Ut 
u A^i A12 ^18 
V ^ij Jg, J^ 

W ilj3 ^23 ilgj 
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those new variables are introduced which would be constants in 
Sir W. R. Hamilton's transformation. 

This remark suggests an extension of the process. If i be a 
function of tf, 0, &c. as well as of ffy <}>, &c. the quantities c^, c^, &c. 
will not be constants. We express this by writing u, v, &c. instead 
of Cj, c,, &c. Suppose we wish to eliminate some of the diflferential 
coefficients, viz. ff, <f}, &c. and to retain the remaining ones, viz. 
^, r/, &c. If we put 

L=^L-u0'-v<f>-&c. 

we may easily show as in the preceding page that 

d dll dL ^ g ^ 

where A is the determinant, called the discriminant, which may he formed from 
the determinant just written down by omitting the first row and the first column. 

This is a general expression for the Hamiltonian function and agrees with that 
which may be dedaced from the result in Art. 21, when all the variables are trans- 
formed by the Hamiltonian process. 

Since the polar reciprocal of the polar reciprocal is the original quadrio, it 
follows that 

^"-d^» ^-57' ^-dS^' 
which are three of the six Hamiltonian equations. 

We may also show geometrically that if the coefficients of T^ be functions of 

dT dT 
any quantity ^, then -^= - -3^» To prove this we notice that it x,yt zhe the 

co-ordinates of a point P, situated on a radius vector OP' of a quadrio 

(0P\« 
Yypf) • The quadrics 

Ti = l and T,=l may be regarded as polar reciprocals of each other with regard 

to a sphere whose radius is fj2 and whose centre is the common centre of the 
two quadrics. Let P be any point on the quadric T^=^l, and let the radius vector 
be produced to Z so that OP . 0Z=2, then the quadric T,=l touches a plane 
drawn through Z perpendicular to OP and Q is the point of contact. Let these 
quadrics be slightly altered in consequence of a variation of 0, so that tiieir 
equations are now Ti-\-dTj^=l and Tm+dT^=l, Let OP and OQ produced cut 
these quadrics respectively in P' and Q, Then 



rx+dn=(^)'. r.+dr.= (gy. 



Now if Z' be a point on OP produced so that OP . OZ'=OP . OZ, the quadric 
7,+c^T, wUl touch the plane drawn through Z' in some point q near Q'. The 
point Q' will therefore lie very nearly in the tangent phuie, so that by similar 
triangles 

OQOZ _ OF 

0Q:~'0Z'~ OP' 

Since each of these ratios is indefinitely nearly equal to unity, it follows that 
dTi:=-dTy 

If we put L=Ti+U and H=sTa- U, Lagrange's equations may be written in 

the forms w' = j^ , »' = jt » ^'=zrr» Hence we have 
d$ dip d}// 

, dH , dH , dH 

-~=^» -^-d0' ""^d^^ 

which are the remaining three of the Hamiltonian equations. 
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We have thus as many equations of the Lagranglan form as 
there are variables f , % &c. Also since ^— js , &c. we have by 

differentiation 

dL' (dL \dff „^. „ 

■dn^yw-Vdii-^^^'^-^' 

with similar equations for 0', &c. By Lagrange's equations we 
obtain 

dL , o 

Thus we have as many sets of equations of the Hamiltonian 
form as there are variables d^ ^, &c.] 

21. We may effect this elimination once for all and find a 
definite expression for L\ 

Let the kinetic energy be /^^^)''4^-'^\ 

T = Tea -jr + T$^ ff<l> + . . . ( — ( \,' ' ' 

Then the integrals used will be r^^-^""'^ ' 

T^ff ^tT^^ + ... = Ci— T^l^ — T^r{ -- ?^ — ^ 

&c. = &c. 

For the sake of brevity let us call the right-hand members of 
these equations Ci — X, Cj— F, &c. Since f is a homogeneous 
function^ we have 

. + K(c. +Z) + 4f («,+ F) + &c. 

If we substitute in the second line the values of ff, <f>', &c. found 
by solving the integrals just written down, we have 




2A 



Cj-Z c^-Y 

Ci"-- ^ 2$$ T$^ 



• • « 
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where A is the discriminant of the terms in T, which contain only 
ffy ^', &c., and may be derived from the determinant just written 
down by omitting the first row and the first column. 

We may expand this determinant and write it in the form 



+ 



2A 










"« 

T, 



$^ 



• •• 



' • • • • I 



2A 



0X7 
X Tit T$^ 



1 
A 







X 

Tti 



Y 



• • • • • 



where X, F, &c. stand for 

&c. = &c. 

The first of these three determinants will contain only the 
constants Cj, c^ &c., and the co-ordinates f, 77, &c. The second will 
not contain Cj, Ci, &c. but will be a quadratic function of f', 77', &c. 
The last determinant will contain terms of the form ^', 17' with 
variable coeflBcients which may also be functions of c^, c,... 

Since f f, 1777', &c. are all small quantities, it is clear that this 
expression for L when expanded will take a form precisely similar 
to that given in Art. 2, only that we have fewer variables to deal 
with. 

22. As an example, let us consider the following problem. 

A body has a point O which is in one of the principal axes at 
the centre of gravity G fixed in space. The body is in steady 
motion rotating with angular velocity n about OG which is verticals 
Find the conditions that the motion may be stable. 

Let OA, OB, 0(7 be the principal axes at and let OG co- 
incide with the vertical OZ in steady motion. Let f, 77 be the 
direction cosinea of the vertical OZ referred to OA, OB. Let 
G>i, G>j, G>3 be the angular velocities about the principal axes at 0, 
Then to the first order 

ft), = 77' + ft>3^] 

o), = - f + 0)377. 
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Let be the angle ZOO, yjt the angle the plane ZOG makes 
with a plane ZOX fixed in space and <f> the angle it makes with 
the plane AOC fixed in the body. Then 

f = — sin^cos0) 
rj=s sin ^ sin ^j ' 

©3 = 0' + '^' cos $ 

=f+(x-*')(i-D. 

putting j^ = + '^. We easily find 

If then J, B, Che the principal moments of inertia at 0, the 
Lagrangian function is 

+^(v+x'fr+|(-r+x'9)' 

where if is the mass of the body, and h » 00^ 

Since ;^ is absent from the equation we have the integral 

dL_ 

which gives 

5^' = n + terms of second order. 

Hence 

L' = L-Cn)l 

= const+^i," + |r 

+ {( J - C) n* + %;!} I + {(S - (7) n' + %M I 

Using this as the Lagrangian function we easily find 



{A - C)r? + Mgh - Bm\ {A + B-C)nm 

-(A + B-C) nm, {B -C)n' + Mgh - Am' 

B. A. 



= 0. 
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The roots of this equation to find m must for stability be of 

the form +/3\/ — 1. Putting m* = — X* we have a quadratic to 
find X*. The roots of this quadratic must be real and positive. 

If -4 = -B, as in the case of a top spinning with its axis vertical, 
we have 

The motion is stable or unstable according as C^n^ is greater or 
less than ^AMgh, If CV = ^AMgh, the equation has equal roots 
and as the first minors are not zero the motion is unstable. 

23. [As another example of the use of the modified Lagrangian 
function, let us consider a case discussed by Prof. Ball in the 
Notices of the Royal Astronomical Society for March, 1877. In a 
problem in Physical Astronomy, we want the relative co-ordinates 
of the system, while its absolute motion in space does not concern 
us. Lagrange's equations involve both the relative and absolute 
co-ordinates, and are therefore not particularly well adapted for 
such problems. By using the modified Lagrangian function, we 
may eliminate the absolute co-ordinates. 

Let the system have n co-ordinates, let us choose as three of 
them the co-ordinates of the centre of gravity of the whole system, 
viz. 0, (f), yfr. There will then remain n — 3 co-ordinates which are 
independent of these. Let T be the kinetic energy of the system 
relative to its centre of gravity, V the potential energy, M the 
whole mass. Then the Lagrangian function is 

In problems in Physical Astronomy the potential energy is a 
function only of the relative positions of the bodies, and is there- 
fore independent of 0, <f>, yp- and their differential coefficients. We 
have therefore 

dL dL dL _ 

dff^^'' df"^^ d^'^^'' 

Hence the modified Lagrangian function is 

L' = T — F— a constant. 

It is independently clear that we might take this as the 
Lagrangian function, for the first three terms of L do not enter into 
any one of the Lagrangian equations, except the three formed by 
differentiating with regard to ff, </>', '^^'^ 

The function T is made up of two parts, (1) the kinetic energies 
of the rotations .of the bodies about their centres of gravity, which 



IV.] LAGRAJ^GIAN FUNCTION. * 67 

we may call T/, and (2) the relative kinetic energies of the several 
bodies, each collected at its centre of gravity, which we may call 
T'. Let m„ m^, &c. be these masses; x., a?,, &c. the abscissae of 
th'eir centred, of Vavity referred to the J^enlre of gravity of the 
whole as origin. Then, accents denoting diflferential coeflBcients 
with regard to the time, we have 

m^x^ + m^x^ + &c. = 0. 
Let us square this and write 

2xX = ^'" + a?.'" - « - <)'. 

If we examine the coeflBcient of any power as a?/" we see that it 
is 

m^ + m^{m^ + WI3 + &c.) = m^m. 

Hence the square becomes 

2m %moc^ — %m^m^{xl — a?/)" = 0. 

Similar expressions hold for the y and z co-ordinates. Hence 
on the whole we see that the relative kinetic energies of the several 
bodies collected at their respective centres of gravity is 



22m ' 

where v is the relative velocity of the centres of gravity of the 
masses m^, m,. If we express this in any kind of co-ordinates, we 
may use the Lagrangian function i' to find the relative motion. 

The expression for T^ agrees with that given by Prof. Ball, but 
his demonstration is quite dififerent. Prof. Cayley has given 
another demonstration in the same number of the Astronomical 
Notices. 

The Lagrangian function thus found may be still further 
"modified." To avoid symbols of summation, let us consider the 
case of three particles moving in one plane under their mutual 
attractions. Let the separate masses be m^, m^ m^, and let fi be 
their sum. Referring the system to m^ as a central mass, let the 
distances of m^ w, from w^ be respectively r and p, and let the 
opposite side of the triangle be R. Let the interior angle be- 
tween r and p be and the exterior angle between r and B be 
X' Let be the angle r makes with some fixed straight line in 
space. We easily find 

where A = mjmy + m{(n^ + mmJS^y 

C = ^r- XmmV, 

Z/Jb 

5—2 
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and V is the relative velocity of the masses m, w! calculated on 
the supposition that mi is fixed, and that the straight line r has 
no rotation round m^. Thus ul, B, C are all functions of r, p, ^ 
and their dififerential coeflBcients with regard to the time. 

If we only want the changes in the form and magnitude of 
the triangle joining the three particles, we may eliminate ff by 
means of the equation 

Aff ^-B^c^. 

We then find as our modified Lagrangian function 

which contains only the three co-ordinates r, p and 0]. 

24. When the geometrical equations contain differential 
coeflBcients of the co-ordinates f , 77, f, &c. of the system with regard 
to the time, we cannot express the co-ordinates a?, y, z of any 
element of a body in terms of |, 17, ^, &c. by means of equations of 
the form 

y =A(.^y Vf 5 &c., t), 

It follows, as is pointed out in our books on Rigid Dynamics, 
that Lagrange's equations cannot be employed in the form 

dt d^ di d^' 
In many of the most interesting problems in Rigid Dynamics, 
it so happens that the geometrical equations do contain —^ 

-^ , &c. For example, let a sphere be set rotating about a vertical 

diameter and be on the summit of a perfectly rough surface of 
any form. If a small disturbance be now given to it, the sphere 
may roll round and round the summit. During this motion the 
velocity of the point of contact is zero, and our mode of re- 
presenting this analytically in terms of the co-ordinates will give us 
two equations of the form 

To include such cases the equations of motion must be modi- 
fied. If L be the difference between the kinetic and potential 
energies, all the Lagrangian equations may be written in the form 

fd dL dL\ 5.^ . /d dL dL\ rs , « ^ 
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TV.] 

where S^, Srj, &c. are any small arbitrary displacements consistent 
with the geometrical equations. But if these geometrical equa- 
tions be given in the form 

&c. = 0- 



these arbitrary displacements must satisfy 

&c 



. = 01 

.=or 



If this were not the case, the geometrical displacement of the 
body given in applying Virtual Velocities would not be such as to 
cause the unknown frictional forces, &c. to disappear. Using the 
principle of Indeterminate Multipliers, we get 

d dL dL ^ xv TT /v 



* » 



d dL dL ^ j^ TT ^ 

dtdi/'d^'^^^^^^^^-^--'' 

&c. = 
These joined to the geometrical equations 

&c.=or 

are suflBcient to determine the unknown co-ordinates f, rj, &c. and 
the multipliers \ fi, &c. 

It will be more convenient to write these equations in the form 
d_dL_dL dG dH 
dtd^' di'^^d^'^^W'^'"^ 
ddL dL^^dO dH ^ _^ 

&c. = 0, 
the geometrical equations being 

G = 0, ir= 0, &c. 

It is of course obvious that these indeterminate coefficients 
\, /Lfc, &c. are merely the frictions or other resistances introduced 
into the equations in a convenient form. 

25. In order to apply these equations to the oscillations of a 
system about a state of steady motion, it will be convenient to 



70 



INDETERMINATE MULTIPLIERS. 



[chap. 



change the co-ordinates f , 77, &c. into others tf, ^, &c. which vanish in 
the steady motion. Let L be thus expanded in powers of 6, <f}, &c. 
as explained in Art. 1, and let P and Q have the meaning given 
to them in Art. 5. 

Let us then put 

|' = a + ^, V = /8 + <^', &c. 

where a, 13, \, fi^, &c. are the values of ^, rf, \ fi, &c. in steady 
motion. The geometrical equations will then take the form 



&c 



.=01 



and the equations connecting Sd, B^, &c. will be 



OJSd + G,Bif> + «fec. = 

&c 



.=0) 

•=or 



In these equations G^, G^j &c. are functions of 0, 0, &c. Let 
a square bracket indicate that the value of the inscribed quantity 
in steady motion is to be taken. Thus [GJ means the value of G^^ 
when 0, 0, &a have all been put zero. 

The equations of steady motion may then, exactly as in Art. 3, 
be written 



-[y+^o[5|] 



+ /^c 



dif/^ 



+ &c. = 
&c. = 0. 



From these the relations which exist between a, /3, &c., Xo> A^o> ^^* 
may be found. 



The equations of the oscillatory motion may be written 
with similar equations for 0, -^j &c. 



=0, 



26. The final determinant written for two variables 0, <f>, and 
two geometrical equations G and H in the notation of Art. 19, 
will be 
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! B,y - A,, + E,, 

'^^'Iddy'^'ldtij 


B^m?-A^ + E^ 
+ {G„-CJm 


-dG' 'dH 
\_dff\ \_de'\ 

'dG' 'dH 
L#'J UfJ 





B„m' - ^„ + ^1, 

^^'lddy'''>ldd] 


B„m* -A„ + E„ 


[dOi ^ [dcn 
[dej^ IdO'i *" 




Idej^ldffi'^ 


"dfl" .'dH' ^ 



= 0. 



It will be noticed how very much this determinant is simplified 
if the values of \, ^ in steady motion are zero. 

27. Let us apply these equations to the solution of the follow- 
ing problem. 

A heavy sphere rotating about a vertical diameter rests in equi- 
librium on the summit of a perfectly rough surface and being 
slightly disturbed makes small oscillations, find the periods. 

As the sphere moves about, its centre always lies on a surface 
which may be called parallel to the given surface. Let the high- 
est point of this surface be taken as the origin and let the axes of 
00 and y be the tangents to its lines of curvature at 0, so that, the 
equation to the surface in the neighbourhood of is 



\Pi PJ 



Let P be the centre of the sphere, PC that diameter which is 
vertical when the sphere is in equilibrium on the summit. Let 
PA J PB be two other diameters forming with PC a system of 
rectangular axes fixed in the sphere. Let the inclination of PC 
to the axis of Z, which is vertical, be 6, and let the vertical plane 
through PC make with the plane scz an angle '^, and with the 
plane CPA an angle ^. The vis viva 2T of the sphere will then 
be given by 

r= i \a)'' + y") + i Z?* {(f + f cos ey + 0" + sin^e^'^'}. 
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Let sin 6 cos -^ = f , sin ^ sin -^ = iy, then we have to the 
necessary degree of approximation 



tf'» + sin»e^'» = ^'^ + 77'»r 



Also let A 4- -^ = ;^. These transformations of co-ordinates are 
all permissible, because they do not involve any differential coeflS- 
cients with regard to the time. We thus find if ii be the differ- 
ence between the kinetic and potential energies 

x=i(«;''+/»)+f{x''-x'(^?'-i?r)+r+';'*}+f^+^. 

If 0)^., ft>y, G), are the angular velocities of the sphere about 
parallels to the axes, the geometrical conditions are 



x' — a {tOy — ft), — j = 
y' +a\(o^-w^^ =0 



where a is the radius of the sphere. These equations by weU- 
known rules reduce to 

1- ^' sin -^ sin tf + ^ cos -^ — ['^' + <^' cos 5) — = 

-^ -d) cos'drsme + ff smylr + Nr + 6' cos 5) - = 
a ^ ^ ^ Pi 

expressing these in terms of our new co-ordinates we have 



Pi 



a 



The position of the system has now been expressed in terms 
of such co-ordinates, that the coeflBcients in the governing func- 
tions L, O, H are all constant. See Art. 2. 

The steady motion is given by x, y, f, 77 all zero, and x — w« 
To find \ fi^ we may use the equations of steady motion 

-[f]--[i]-.[f]=«' 

where q stands for any one of the co-ordinates. Taking j = a? and 
gr = y, we see that Xq = and ^^ = 0. 

To find the oscillation we may substitute in the determinant 
and thus form the equation which gives the periods. As most of 
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the constituents of the determinants are zero, it will be more con- 
venient to form each equation directly from the standard formula 



dq 



^--g-.[f]-[f]=°. 



where q stands for any one of the co-ordinates. Taking q in turn 
to be X, y, x, f, rj we find 

Pi a 
^ ^P> « 

^(r+xV)+\=o 

^ iv" -XVI + 1^, = 

Putting x'=^ these with the two geometrical equations are all 
linear and ready for elimination. 

Eliminating \, fi^ we have 



k 



Afp. 



Substituting for f, rj from the two geometrical equations, we 
have 



Pi <^P. 



a' 



a />2 » Pi 



► • 



r-« 



To solve these put x = Xcos (pt +q)y= Fsin (pt + q) so that 
/? is the quantity required. We obviously have 

\r ^ a' + i? pj [^ ^ a' + il' pj'' (a' + J^r pj>/' 
which is a quadratic to find p\ 

If p^y p^ have opposite signs the roots cannot be real, and the 
steady motion naust be unstable. If p^ p^ are both positive, so 
that the sphere is on the summit, the motion is stable only if 



n > 




CHAPTER V, 



Certain euhaidw/ry determinants are formed from the dynamical de- 
terminant, and it it shown that there m,uel be at least as many 
roots indicating stability as there are variations of si^ lost in 
these subsidiary determinants, and must exceed the number lost by 
an even number. Arts. 1 — 5, and 9. 

2ftt8 M equivalent h 

with simila/r limitations. 

Effect of equcU roots on this test of Ae ttabUity of the system. Arts. 
10, 11. 

Example, Art. 12. 

1. Id order to test the nature of the roots of the determi- 
nantal equation, let U3 apply a method analogous to that by which 
Dr Salmon in his Higher Algebra proves the reality of the roots 
of the equation which occurs in the determinatioa of the secular 

inequalities of the planets. 

2. Let A he the determinant which forms the left-hand side 
of the fundamental equation, let A„ he the determinant formed 
by omitting the ^* row and 5* column. Let A, be the determi- 
nant formed by omitting the first r rows and r columns. Thus 
A, = A„. We then have by a known theorem in determinants 

AA, = A„A„-A„A„. 

It has already been noticed that if we change nt into — m, the 
determinant A is changed into another determinant whose columns 
and rows are the same as the rows and columns of the first de- 
terminant. It easily follows that the minor —A,, is chained into 
the minor — A„ by changing m into — m. Hence if 

A„ = <^K) + m>>-(m')) 
then A^ = (m') - m-j- {n^)) ' 
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Hence the product Ajj,Agj is necessarily positive for all negative 
values oT m*. 

• • • - 

It also follows that if A * y^isheg for any negative value of n^ 
then A„ also vanishes for thfi iaiS^ Value of w*. 

3. When the determiiiant/A; vanishes, we have 

so that A and A, must have opposite signs, or one of them must 
be zero. Consider then the series of determinants 

A, Aj, A„ A„ 

each one being formed from the preceding by erasing the first 
row and the first column. We thus have a series of functions 
of m* whose degrees regularly diminish from 1ihe ifi^ to the first. 
As we may suppose the determinant A to have a row and a column 
of zeros added on at the bottom and right-hand side, but with any 
positive constant in the right-hand bottom comer, we may add 
to this series of determinants any positive constant. We have just 
proved that if any determinant of this series vanish for a nega- 
tive value of m*, the two determinants on each side have opposite 
signs. The case in which two successive determinants vanish for 
the same value of m' will be considered afterwards. 

We may then use these determinants in a manner somewhat 
similar to that in which we use Sturm's functions, provided no two 
successive functions vanish for the same negative value of m*. 
No variation of sign can be lost as we pass from m^ = — oo to 
m' = except by the vanishing of the determinant A at the head 
of the series. And when a variation of sign is lost, it will be 
regained again at the next root, unless a root of the determinant A^ 
separates the two roots of the determinant A. If therefore in this 
passage from — oo to zero, as many variations of sign are lost 
as is indicated by the highest power of m*, the valves of m* fownd 
from the determinantal equation must be all real and negative. 
It will also follow that the roots of each of the series of determi- 
nants are all real and negative, and that the roots of each separate 
the roots of the determinant next above it. If, however, the 
proper number of variations of signs be not lost in the passage 
from m' = — 00 to m* = 0, it does not follow that the values of m' 
are not real and negative. 

4. If the proper number of variations of sign has not been 
lost in this passage from m' = — x to m' = 0, this proposition 
does not leave us without information as to the nature of the roots. 
We infer that the number of real negative values of m' is equal 
to or exceeds the number of variations of sign lost by an even 
number and unless tlie number of variations be even the system 
is unstable. 
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If we do not object to the labour of expanding the determi- 
nants, we might extend this theorem to determine the positions 
of the negative values of m^ as well as their number. The number 
of real negative values of m' between m* = — a and m* = — jS is 
equal to, or exceeds by an even number, the number of variations 
of sign lost in the series of determinants. In this form the 
theorem resembles Fourier's theorem in the Theory of Equations. 

5. The converse of this proposition has not been proved. 
If in the passage from m* = — oo to m' = 0, the number of varia- 
tions of signs is unaltered, it is not true that the values of m' 
cannot be real and negative. Thus in the simple case 






0, 



where a^^, a^,, a^ are all positive, no variations are lost, yet if 
^12 > ^^11 "^ ^^M ^^^ values of m' are real and negative. And if 
^12==^'^ + ^^ ^^ roots are equal and negative. It will be 
noticed that the minors in this last case are not zero. 

6. In order to discover the meaning of these losses or gains 
of changes of sign, it will be convenient to make such changes of 
the co-ordinates as will simplify the dynamical determinant as 
much as possible. Let us write 

If we now change our co-ordinates by writing for ^, ^, &c. 
linear expressions of some new co-ordinates, we know that we can 
clear this expression of all the terms containing the products. 
We know also that this can be done in an infinite number of 
ways. We thus have 



'o ~ ^11 "2" "^ ^22 "2" + • ' * 
where the symbols 0^, ^^, &c. represent the new co-ordinates. 

Again, let us consider the expression 

TF. = 5,^ + J5,,^f + 

the coefficients are here the values of P, Q, &c. in the general 
expression for T, Art. 1, Chap. iv. when 9, <^, &c. are all put zero. 
But since T is necessarily positive for all values off', rj'...^, rj,... 
it follows that W^ is positive for all values of ff, <f> ... Hence 
by a well-known theorem, we may by a real linear transform- 
ation of the variables clear the expression W^^ also of the terms 
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containing the products 6\ ^', &c. and can make the coefficients 
of the squares any positive constants we may please. We thus 
have 

Tf =^ + ^^. 

^U 2 2 



It may be shown that we cannot in general clear the expres- 
sion for T of all the terms containing 6^^ a if), &c. unless G^^— 0,^=0, 
&c. by substituting for 0, ^, &c., any linear functions of other 
variables. As this is only a negative result of which no further 
use will here be made, it is unnecessary to supply the demon- 
stration. 

7. Using these simplifications the determinant A will now 
take the simpler form, 



A=: 



«12^> 


a„w. 


^' + CL^, 


a^m, 




m' + a33, 
> 



= 0, 



> 

where a^ = (7^ - C^^ &c. 

If we form from A the series of subsidiary determinants 

A, Ai, A,..., 

terminating with any positive constant, we see that when w' = — oo, 
these subsidiary determinants are alternately positive and negative, 
and when m^ = 0, they become 

which are all positive if 0^, a^... are all positive. Hence if 
ctu> ^sw« are all positive, the proper number of changes of sign has 
been lost, and therefore the roots of the dynamical equation A = 
satisfy the condition of stability. 

If Oil, a^, &c. are all positive, we see that V^ is a minimum for 
all variations of 9^, ^j, &c. and therefore for all variations of the 
original co-ordinates. If a^^, a^^, &c. are not all positive, there will 
be as many variations of sign lost as there are positive quantities 
in the series Oj^, a^, &c. In this case V^ is a minimum for some 
variations of u, ^, &c. and not for others. It is shown in the 
appendix to Williamson's Differential Calculus, that n independent 
conditions are necessary that a quadratic expression of n variables 
should be always positive. These are given in the form of deter- 
minants and may be briefly summed up, in the statement that 
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= 0, 



is always positive if the roots of 








-a.jj + A,, -^12 > 


-^18» 


<fec. 




■^12> -^89 "^ ^» 


^28 > 


&c. 




^18 » -^23 > 


As + \ 


&c. 




&C. &C. 


&C. 


&c. 


are all real and 


negative. 







8. We may now put the proposition of Art 3 into another 
form. Let L be the general expression for the excess of the kinetic 
energy over the potential energy of a dynamical system in terms of 
its n co-ordinates x, y, &c. Let this system be moving in steady 

motion with constant values ^/ jt > jf > &C' Then if h be a rnaxi- 

dx dv 
mum for all variations of x, y, &c. keeping -ri * -A> &c. unchanged, 

then that steady motion is stable. 

If however only r of the n conditions necessary to mxike L a 
ma^mum be satisfiedy then the number of roots of the dynamical 
equation which satisfy the conditions of stability is equal to r or 
exceeds r by an even number. There cannot be stability unless n — r 
is an even number. 

If the system be oscillating about a position of equilibrium, 

dx du 

-jif-Ty &c. are all zero, and this leads at once to the condition, 

at at 

that the equilibrium is stable if the potential energy is a minimum. 

9. In this reasoning, we have for convenience excepted the 
case in which two successive determinants in the series A, A,, A,. . . 
vanish for the same value of w*. But this exception is of no real 
importance, for we may change these determinants into others 
whose constituents are very slightly different from those of the 
given determinants but which are such that no successive two of 
the series have a common root. In the limit, therefore, when 
these arbitrary changes of the constituents are indefinitely small, 
the roots of the series of determinants will still be real under the 
same circumstances as before, and the roots of each will separate, 
or coincide with, the roots of the next above it in the series. 

To show that these changes are possible, let us consider the 
row of determinants beginning at the last. The determinant A„ is 
a positive constant, the next A,^^ is m* + a^^. Proceeding thus, 
suppose we arrive at two determinants which we may call A^ and A 
which have a common root. If we now change the constituents 
Ou, ttjj, a^g, &c. into a^ + ha^, a^^ 4 ha^^, &c. we do not alter A^, but, 
except for the root w* = 0, we do alter A in an arbitrary manner. 
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When for example a,g is altered, we alter a constituent both in the 
first row and in the first column. Since 

AA' = AA -A A , 

where A' is the determinant formed from A by omitting the first 
and third rows and columns, we see that when A and A^ both 
vanish, the product AigAji and therefore by Art. 2 both Ajg and 
Aji must vanish. The determinant A is therefore altered by 
— Ag (Sajg)' which does not vanish, since A, is by hypothesis finite 

for the particular value of rr? under consideration. 

♦ 

If any determinant of the series vanishes when rn? = 0, it is 
clear that one of the quantities a^, a^.,, must be zero. If we 
replace this by any small positive quantity, the argument will 
apply as before. 

10. It is important to consider the effect of equal roots on the 
test of stability given in Art. 8 of this chapter. 

In this case we know that the roots of the minor A^ separate 
the roots of A. K therefore A have two negative equal roots, it is 
clear that Au must have one of them. In the same way Ajj, A^g, 
&c. must also all vanish for this value of m. Since 

AA' = A A -A A . 

it follows as in Art. 2 if A and A^ both vanish, that A^^ and A^^ 
also vanish. Hence all the first minors of the determinant A 
vanish. This is the case considered in Art. 5 of Chap. I. The 
equal roots instead of introducing into the expressions for the co- 
ordinates terms which contain ^ as a factor merely render two of 
the coefficients, instead of one, indeterminate. 

In the same way if the equation A = is satisfied by three 
values of m^ equal to the same negative quantity, the equation 
A^ = must have two of them, and its principal minor must have 
one of them. Reasoning as before we see that all the second 
minors of A must be zero. This is again the test that there should 
be no terms which contain.^ as a factor. 

The presence therefore of equal roots does not in the theorem 
of Art. 8 affect the stability of the motion. 

When a system is disturbed from a position of equilibrium 
whether stable or unstable, the roots of the fundamental determi- 
nant are separated by its minors in the manner described in 
Art. 8 of this Chapter. By what has just been proved, we see 
that if the fundamental determinant have equal roots, whether 
positive or negative, these do not introduce into the integrals 
terms which contain < as a factor. 
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11. The proposition in the last article may be made more 
general. If the fundamental determinant be reduced to the form 
indicated in Art. 7, we shall show that if A vanish for two equal 
negative values of m' which are numerically greater than the 
greatest negative quantity in the series a^^, a^, a^, &c., then 
these equal roots will not introduce any terrifis into the solution 
with ^ as a factor. If a^^, a^y &c. are all positive, this reduces to 
the proposition proved in the last article. 

Following the same dotation as before, we have 

AA' = A A -A A . 

If neither A^^ nor A„ are zero, they must have the same sign 
when A vanishes for a negative value of m\ For their product 
is equal to AjjA,^ which has been proved in Art. 2 to be positive. 
Hence all the leading first minors, viz. Aj., A^, &c. must have the 
same sign for any negative value of m^ which makes A vanish. 

By differentiation we have 

^ = 2w2A„ + a«Ai2+--«i2^2i + 2mA,, + ...-&c. 

But we have also 

A = (w" + a^^) Ajj + a^m A„ + ... 
A = -a„m A„+ (w' + ajj Ajj + ..^ 
&c. = &c. 

Hence if w be the highest power of m occurring in A, we have 

/7A 
»» ^ = «^ + (»»*-««) Au + ("*'-« J ^ - &c- 

/7A 
If then A and -j— both vanish for any negative value of w* 

greater than the greatest negative quantity in the series a^^, a^^, 
&c., we have the sum of a number of quantities all of the same 
sign equal to zero. This requires that each should be zero. We 
have therefore \^==0, A^=0, &c. The rest of the proof is the 
same as before. 

/7A 
By differentiating the expression for -j— and substituting for 

J A /7A 

-^" , --T-^' , &c. their values in terms of their first leading minors, 

we may extend this proposition to the case in which the funda- 
mental determinant has three equal roots, and so on. 
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12. A sphere is suspended by a string OA from a fixed paint 
0, and is set rotating about a vertical diameter which is in the 
same straight line as the string with an angular velocity n. A 
small disturbance is given, determine if the steady motion is 
stable. 

Let be the origin, and let the axis of z be vertically down- 
wards, let Ix, ly, I be the co-ordinates of -4, the point at which the 
string is attached. Let G be the centre, and let af, at), a be the 
co-ordinates of C relative to A. Then, exactly as in Chap. iv. Art. 
22, the Lagrangian function may be shown to be 

the mass being taken unity. 

Putting x' = ^y ^f Vf ^\ y' ^ zero, we see that' Z is a maxi- 
mum when Xy y, f, 17 are zero, the steady motion is therefore stable 
for all values of w» 

If we put V — \ofy and m' = — X', so that x, y, &c. are all repre- 
sented by terms of the form 2-4 cos (\^ -f- a), we may, by the 
methods of the last chapter, prove 



i-T 



(^■-D^-i-^-feM^- 



This equation, whatever may be the sign of w, has two positive 
and two negative roots. All four give stable oscillations. 



R.A. 



CHAPTER VI. 



If ilie energy of the system he a maadmum or minimum under certain 
conditions, the motion whether steady or not is stable. Arts. 1 — 3. 

When the motion is steady, it wiU he also stable if a certain funcUon of 
the co-ordinates called Y + W w a minimum. Art. 4. 

If there he only one co-ordinate which enters into the Lagrangian function, 
except as a differential coefficient, this condition is necessary and 
sufficient. Arts. 5, 6. 

Additional conditions when there are tvx> co-ordinates^ Art. 8. 

1. When a system is oscillating about a position of equilibrium, 
it is well known that we may determine the stability or insta- 
bility of the equilibrium by what we may call the " energy crite- 
rion." This criterion may also be sometimes used when the 
system is oscillating about a state of steady motion. 

Let E be the sum of the kinetic and potential energies of the 
system. Then throughout any motion of the system we have 

E=^h, 

where A is a constant depending on the initial conditions. If 6, (f>, 
&c. are the co-ordinates of the system, J? is a known function of 
6, , if), <j>',&c. Suppose that some of the other first integrals of the 
equations of motion are known. Let these be 

F, {0, ff, <l>, ^\ &c.) = G, 
F, ifi, ff, </>, 4>\ &c.) = C, 
&c. 
the time t being absent. 

For the purposes of this proposition let us suppose 6, ff, ^, ^', 
&c. to be separate variables unconnected with each other except 
by the equations just written down. 

If E be an absolute maximum or an absolute minimum for 
all variations of 0, ff, &c., those corresponding to the given motion 
making E constant, then that motion is stable for all displacements 
which do not alter the constants (7,, CL, &c. 
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If this proposition be not evident, it may be proved by elimi- 
nating as many of the letters as possible. If ^, 0\ &c. be the re- 
maining co-ordinates we have 

Let h be the value of E in the given motion, and let the 
system be started in some slightly different manner so that 

If J? be a maximum along the given motion, then any change 
whatever in 6, 0\ &c. decreases E, Hence 0, d\ &c. cannot deviate 
so much from their values along the given motion that the change 
in E becomes greater than Bh. 

2. Let us apply this principle to a system of bodies which 
moves in steady motion with some co-ordinates 0, <f), &c* such that 
their differential coefficients 0', <f>', &c. are constant, and the re- 
maining co-ordinates ^, rj, &c. themselves constant. Let us further 
suppose that the energy is a function of 0'y <f>\ &c., but not of 
0f <f>. By Lagrange's Equations we have the integrals 

dT_ ^ dT __ ^ o 
dff" ^' d<i>'' ^' 

It is clear that the system can describe any one of a number 
of steady motions, which we have already called parallel motions, 
and which are determined by 

ff =p, <!>' = q, &c. 

f = af, 97 = /3;&c. 

where p, q, &c. a, ^, &a are constants which satisfy all Lagrange's 
equations. 

We have thus as lioany relations between these constants as 
there are co-ordinates f, rj, &c. Let the system be started with 
any initial conditions we please, then the constants C^, (7,, ... are 
given. These being known we have as many relations between 
the constants of steady motion as there are co-ordinates 0, if>, &c. 
The steady motion is therefore determined. • If the energy of this 
initial motion is nearly equal to that of this steady motion, and 
if it be a maximum or minimum as explained above, then the 
system will never deviate far from its corresponding position in 
the steady motion, and this steady motion may be called stable. 

3. Example. A top is set spinning on its point on a perfectly 
rough horizontal ground, with its axis inolined to the vertical, find 
the condition of stability. 

Let be the inclination of the axis OG of the top to the 
vertical OZ; i/r the inclination of the plane ZOC to a vertical 

G— 2 
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plane fixed in space^ and ^ the inclination to a plane through OG 
fixed in the body. Let be the apex, the centre of gravity 
which lies in OG,h=^ OG. Let A, A, G be the principal mo- 
ments of inertia at 0, and Jf the mass of the top. We have then 

E= ^ (f + ^^' cos ^» + ^ {ff* + sin'^^'«) + Mgh cos ^. 

By Lagrange's equations we have the integrals 

^' +'^'cos^ = n, 
Cn COS0 + A sin*^^' = m, 

where n and m are two constants, the former representing the 
angular velocity of the top about its axis, and the latter the 
angular momentum about the vertical. If we now eliminate 
<l> and y^' we find that J? is a minimum when ^ = a^ if 

Cn' > iMghA cos a, 

which is the result given by other methods. 

4. The theorem of Art. 2 may be put into another form. 
Let the kinetic energy be 

T= 2i«-2-4- T$^0^ + ... 

Then since 6^ 6, &c. are absent from the coefficients^ we have 
the integrals 

Tees' + ?Wf + ... = a,- Te^^ - ^^V - &c. 

3V^^ + r^f + ... = 0, - 2wr - T^V - &c 

&c.=&c. 

For the sake of brevity let us call the right-hand sides of these 
equations G^ — X, C, — Y, &c. Since T is a homogeneous func- 
tion of ff, ^', &c., we have, as in Chap. IV. Art. 21, 

+ ^ff{G, + X) + i<l/{G,+ Y)+... 

If we substitute in the second line the values of ff, ^\ &c. 
found by solving the linear equations just written down, we have 
the determinant 



0, G^+x, c,+ r, &c. 

(7,-X, Tee, ?W> &c. 

c,-r, r^, r^, &c. 

&c. &c. &c. &c. 



1^ 

2A 
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where A is the discriminant of the terms in T which contain 6\ 6\ 
&c. This determinant is unaltered by changing the signs of X, i, 
&c. and is a quadratic function of C^, C,, &c., X, Y, &c. Hence 
the teims C^X, CjF, &c. do not occur. If then we put 



F=- 






c., 


c^.> 


&c. 


Oo 


Te0, 


T^. 


&c. 


G,, 


Te^i 


T^, 


&c. 


&C. 


&c. 


&c. 


&c. 



2A 



we have 



T^B^X + B^^v' + Scc+W, 



where B^^, &c. are independent of G^, (7„ &c. Now T is essentially 
positive for all values of the variables, and therefore for such as 
make C^, C^ &c. all zero. Hence the terms involving f ', 17', &c. 
are together a minimum when f , rf, &c. are all zero. The coeffi- 
cients -Bff, &c. may all be treated' as constants since f, r/, &c. are 
all small quantitiies. 

If V be the potential energy, we have therefore the following 
rule. If W + Y be a minimwni for all variations of f , ^, &a then 
the steady, motion i& certainly stable. It, should be noticed that 
W+ F is a function only of f, 17, &c. the corordinates which are 
constant in the steady motion. 

5. If the energy>he a function of one only of the co-ordinates, 
though the differential coefficients of all the omers enter into its 
value, this condition is suffwient and necessary. 



Let f be this co-ordinate. Then by vis viva we have 

2 

Differentiating we have 



■B»Tr + F+F=fe. 



This equation must be satisfied by the steady motion repre- 

d(W+ F) 
sented by f = a. The second term ^ ,^ must therefore 

of . 

vanish when f = a, so that TF+ Fis a maximum, or minimum. 

To fiud the oscillation let us put f = a + ^, we find 
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where the square bracket implies that f is to be put equal to a 
after diflFerentiation. By the same reasoning as before B^^ is 
necessarily positive, and the motion will be stable or unstable ac- 
cording as (W + Y) is a minirnum or Tnaadmum, 

6. If we refer to Art. 21 of Chap. iv. we see that this function 
Tr+ Fis the value of the modified Lagrangian function L' when 
f, 7j', &c. are all put zero, and the sign of the whole function 
changed. It therefore follows by Chap. v. that when W+ V is a, 
minimum the steady motion is stable. The "energy criterion of 
stability," as far as it applies to steady motion, may therefore be 
deduced from that given in Chap. V. Art. 8. But the mode of 
demonstration adopted in that chapter gives us more information 
as the nature of the motion, while the modes of application to 
examples of the two criteria are quite different. The energy 
criterion may also be sometimes applied to determine the stability 
of a motion which is not steady. 

[The relation between the theorem in this Chapter in which 
E= T+ V is made a minimum to that given in Chapters iv. and v. 
in which i= T— F is made a maximum may be more distinctly 
perceived by the following statement. 

Let a?, y, &c. be the co-ordinates of the system, and let L be 
the Lagrangian function* so that L = T— V, then by Art. 8 of 
Chap. IV. the co-ordinates in steady motion satisfy the equations 

Here L is expressed as a function of x, y, &c. x\ y\ &c. 

Suppose some of the co-ordinates as 6, <f>, &c. are absent from 
the expression for L, so that i is a function of 0\ ^', &c., the re- 
maining co-ordinates, viz. f , 77, &c. and their differential coeflScients. 
Then if we form the modified Lagrangian function as in Art. 21 
of Chap. IV. the equations (1) of steady motion become 

f-». f=»-*« (2)- 

Here, as in the Hamilton ian equations C^, C^, &c. are the 
0, (f), &c. components of momentum, and L^ is expressed as a 
function of C^ , 6\, &c. |, rj, &c. |', 97', &c. 

dT 

* Let M, V, &c. be the rr, y, &c. components of momentum, so that m= — ,, &c. 

ax 

and let H be the Hamiltonian function. Then ff=T+ V and we easily deduce 

from the Hamiltonian equations that, in steady motion, 

—- = 0, , =0, &c. 
dx ay 

Here H is expressed as a function of x, t/, &c. u, v, &c. 
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But . L' = T - C,ff - C,<l>' - &c. - V 

by Enler's theorem on homogeneous functions. In steady motion 
f , Tj, &c. all vanish, hence the equations (2) become 

df "' dv ""'^^ w> 

where W is the value of T when ^, 17', &c. are all put equal tq 
zero. Here TFis expressed as a function of 0^, C,, &c. f, ^, &c. 

It is shown in Chap. v. that if the Lagrangian function ex- 
pressed as required in equations (1) be a maximum the motion 
is stable. It is shown in this Chapter that if the function W+ V 
be a minimum the motion is stable.] 

7. To find the condition of stability when the Lagrangian 
function is a function of two only of the co-ordinates, though the 
differential coefficients of all the others enter into its value. 

Let f, ff be these two co-ordinates, then the modified Lagran- 
gian function as explained in Art. 20 of Chap. iv. will be a func- 
tion of f, 17, f , f/ only. 

Let the steady motion be given by f=a, i; = /9, with the 
corresponding values of the other co-ordinates 0, (f>, &c. Then 05 
and )8 are constants. Let ^ = (i + oo, 7j = ^+y, and let us expand 
the modified Lagrangian function in powers of os, y. Neglecting 
the terms of the first order, as they only give the steady motion, 
let 

Also let -&= Cfl— Cjs. Then the condition of stability is that 
the roots of the following equation should be of the form ^J— 1. 



B^rrj?^A^^ B^rrf — A^-\' Em 
B^nf — A^ — Em, B^m* — A^ 



= 0. 



If 



A' = ^u^a - ^la' 1 

and = A^B„ + A^B^ - 2A^B^ 
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be the two discriminants and the other invariant, this leads at 
once to the conclusion that the motion is stable only when 

(1) A is positive. 



(2) ^- is positive and > 2 VA A'. 

If A^-^ + Ai^ + -4jj |- is a maximum when x and y iore zero, 
the two conditions are obviously satisfied. 

This condition may be otherwise expressed ; if i' be the modi- 
fied Lagrangian function, a steady motion is given by 

f =0. f .0, r=o. ,--0. 

This motion will be stable if for the values of £ i; thus found, 

(d^v cTry ( d^L'd^v d^ud^u ^ d?r d}L \ 

^^^ \d^dr{ d^dv) tdf df)'' "^ dv' d^ ^d^dvd^'dv'] 
is positive and greater than 

^(d^LdT d'L']i ( d'L d^ii d^r \ I 

^\dl^ dv* d^dv] '\d^^ dri\ d^'dri] * 

8. The nature of the motion when thus reduced to depend on 
two co-ordinates may be illustrated by geometrical reasoning. Let 
the position be defined by two co-ordinates x, y which are zero 
along the steady motion. Let these be regarded as the co- 
ordinates of a poiiit P referred to any axes. Ihen the motion of 
P exactly represents that of the system. 

Let us construct the comics 

Then, exactly as in Arts. 11 — 14 of Chap. I v., it will be found 
convenient to transform the co-ordinates by writing 



x=^a^p + \q\ 



If fi be the modulus of transformation we have fi = aj)^ — ajb^. 
It is easy to see by actual substitution, if E' be the difierence of 
the coefficients of i)^ and pq, that 
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If the transformation be from one set of oblique co-ordinates 
to another, let od, to be the angles between the axes. We then 
have 

E' ^ E 

sin (o' sin a> ' 

Transforming the axes to the common conjugate diameters, the 
conies become 



XL a * -^a 



r 9 



•^u 2 ^ « 2 "" 

the signs of a and h being such as to make these conies real. The 
equation to find m becomes 

It is therefore necessary for stability that the conic a should be an 
ellipse as well as the conic K It is also necessary that 

both roots having the same sign and the inequality being nu- 
mericaL 

Let OP, OP ; OQr OQ hQ the common conjugates of the two 
conies, this condition then becomes 

E area of conic 6 OQ OQ 



sino)" ^Va6 OP^ OF' 

When the system describes an oscillation with o?ie period, {,e, 
an harmonic oscillation, the path of the representative particle is 
easily seen to be 

(BJW - ^«")/ + (^m"^" - ^« ') 3* = constant. 

The harmonic paths are therefore ellipses. It also appears that 
the two ellipses which represent the two harmonic vibrations and 
the two ellipses a and h have, aU four, a common set of conjugate 
diameters. 



CHAPTER VII. 



Any small term of a high order, if its period is nearly the same ns thai uf 
an osciMation of the systefm, may prqdiLce important effectg an the 
magnitude of the osciMation, Art 1. 

Origin of such tenmsy vnth an example. Arts, 2 — 3. 

Supposing the roots of the determinantaX equation to satisfy the conditions 
of stability to a first approximation, yet if a comm^nsv^ahle relation 
hold between these roots it is necessary to examine certain terras of the 
higher orders to determine whelher they wUl idtimately destroy the 
stability of the system^ Art. 4. 

If a certain relaUon hold amwng the coejfflcients of these terms, they will 
not affect the stability of the system, but only slightly alter the periods 
of oscillation. Arts. 5 — 7. 

Examples, the first taken from Lagrange's method of fmding the oscilkh 
tions about a position of equilibrium. Arts. 8—9. 

If the coefficients of the equations of motion should not be strictly conr 
stant, but only nearly so, the stability vnll not be affected, unless the 
reciprocals of their periods have commensurable relations with the 
reciprocals of the periods of oscillation pf the sy stein. Art, 10. 

1. If we understand that a motion is called stable when any- 
small disturbance does not cause the system to deviate far from its 
undisturbed motion, it is clear that we cannot be certain of the 
stability without examining the terms of the second order. It is 
possible that some of these may have their periods so timed that 
their effects may accumulate until the motion is changed. 

Keturning to the equations of Art. 3, Chap. I. we shall have on 
the right hand, instead of zero, a series of small terms of orders 
higher than the first. To find a second approximation, we substi- 
tute the values of x, y, &c. given at the end of Art. 3, in these 
terms. 

They will therefore take the form Ne"** and will produce in 

(T, y, &c. terms of the form -^j-r e***, where N' is of the same order 

at least as the term considered, and/(m) has the same meaning as 
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in Chap. I. These will have to be expressed in trigonometrical 
real forms, but it is unnecessary to exhibit the process, for we see 
at once that no small term or force (whatever it may be called) of 
a high order can aflfect the stability of the motion unless it makes 
f{n) very nearly or exactly equal to zero. In this case its period 
is very nearly or exactly equal to one of the periods of the motion 
given by taking terms of the first order only. 

A remarkable use of this principle was made by Captain Kater 
in his experiments on the magnitude of gravity. It was important 
to determine if the support of his pendulum was perfectly firm. 
He tells us that he had recourse to a delicate and simple instru- 
ment the sensibility of which was so great that had the slightest 
motion taken place in the support it must have been instantly 
detected. The instrument consists of a steel wire the lower part 
of which, inserted in the piece of brass which serves as its support, 
is flattened so as to form a delicate spring. On the wire a small 
weight slides by means of which it may be made to vibrate in the 
same time as the pendulum to which it is to be applied as a test. 
When thus adjusted it is placed in the material to which the 
pendulum is attached, and should this not be perfectly firm its 
motion will be communicated to the wire, which in a little time 
will accompany the pendulum in its vibrations. This ingenious 
contrivance appeared fully adequate to the purpose for which it 
was employed, and aflforded a satisfactory proof of the stability of 
the point of suspension. See FhiL Trans,, 1818. 

2. Since the term Ne^ is obtained by compounding the 
different terms in the values of a?, y, &c. it is clear that 



• • • 



where p, j, &c. are positive integers whose sum is the order of the 
term. It is therefore only when the roots of the dynamical equa-, 
tion/(m) = are such that a linear relation of the form 

pm^ + qm^ + ... = m^ very nearly 

exists between them, that we may expect to find important 
terms among the higher orders. The order of the terms to be 
examined will he p + q + ..,, and unless this be also small, the terms 
will probably remain insignificant. K the root m^ should occur 
twice in / (m) =0 it is clear that the divisor / (n) will be a small 
quantity of the second order, and the term may be said (as in 
the Lunar Theory) to rise two orders. 

3. To take an example, let us suppose a particle to be describ- 
ing an ellipse about a fixed centre of force in one focus. If dis- 
turbed it will describe a slightly different ellipse. Since 

r = a{l — ecos (n^ + € — ft)) + ...}, 

^ = n< + 6 4- 2e sin (w^ + € — ft)) -f ... 
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we see that a slight change in the elements will cause variations 

in r and 6 of the period — , an additional variation in of the 

^ n 

form ^ + Se and an additional variation in r of the form Sa. All 

these variations should by Art. 3 of Chap. L be indicated by 

expressions of the form 

where the values of m are the roots of the equation /C»w) = 0. 
The roots therefore of the equation /(m) = for hr are m = 

and + W — 1, and for 8^ are wi = 0, and ± n V — 1. We there- 
' fore infer that any small disturbing causes of the second order 
whose periods are nearly equal to that of the particle, will cause 
important inequalities in both 8r and W, and (since / (m) = has 
two roots equal to zero) any tenn of long period will rise two 
orders in iO* 

4. If the roots of the subsidiary equation are such that the 
relation 

holds accurately, the solution changes its character. "We have 
now in the value of a; a term of the form te^^^^ Unless the real 
part of m^. is negative, this indicates that the system will depart 
widely from the motion which we took as a first approximation. 
We must therefore modify our first approximation (as m the Lunar 
Theory) by including in it the terms which produced these im- 
portant effects. We may then, enquire how far this modified first 
approximation indicates that the motion is stable or imstable. 
When these terms are included the equations to be solved are in 
general no longer linear, and it is sometimes impossible to find a 
solution sufficiently accurate to serve as a first approximation 
throughout the whole motion. 

5. In some cases, however, the oscillations may still be 
represented by expressions of the form 

&c., 

where the values of n^, n . . . differ but slightly from the roots of 
the equation /(w) = (X Let us investigate the condition that this 
should be true, and also determine whether the changes in the 
values of m^, m,. . . are sufficient to affect the stability. 

Suppose that we have completed our first approximation, and 
find on proceeding to a higher approximation that the terms 
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present themselves on the right-hand side of the first set of 
equations in Art 3, Chap. I. These terms are supposed to have 
arisen from several relations of the -form 

pm^ + qm^ + ... = wi^. 

If these terms can be included in the first approximation by 
writing n,, tij, &c. for w^, m^, &c. we have, by substitution in the 
diflerential equations, certain equations connecting n, Jf, M\ &c., 
whose left-hand sides are the same as those used in Art. 3 with 
111 written for w-, but on the right-hand sides we have instead of 
zero the quantities N^y N(, &c. The test of the success of the 
process is tnat these modifications in the values of m must satisfy 
the same relations as before. 

Now -Wj, J\Y, &c. are all at least of the second order of small 
quantities, hence up to terms of the first order the ratios M^, 2//, 
&c. will be the same as before, so that we may put 

following the same notation as in Art. 3, Chap. i. We also have 
Let Wj = Wi + Swij we find X\^ 

Similarly 

&c. = &c. 

It is evident, by the theorem of determinants alluded to in 
Art. 3, that these are symmetrical expressions. 

6. We may conveniently express these results in the form of a 
rule. 

Suppose we have to a first approximation 

Eliminate from the differential equations all the variables except 
X in the usual manner. This may be done by performing on the 
several equations the operations represented by the minors a, a, 

&c., -J being written for iw. Let the equation thus found be 



/@a. = J.«-.' + P,-«'... 
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Then all these terms can be included in the first approximation, 
provided 

satisfy the relations 

pBm^ + qBm^ + . ,. = Bm^, 

&c. = &c. 
which exist among the roots of the dynamical equation. 

7. The general results we have arrived at may be summed 
up as follows. Though some of the terms of the higher orders 
may affect the magnitude of the oscillation, yet no term will arise 
to affect the stability of the motion unless there be some relations 
between the roots of the dynamical equation of the form 

where p, q, &c. are all integers. Even if such relations occur, the 
lowest order of the term is ^ + ?+ ..., and if this be considerable 
the term will not produce any important effects until a con- 
siderable time has elapsed. If a certain relation, just found, hold 
among these terms, their only effect is slightly to modify the 
periods of oscillation, without altering the type of motion. 

8. As an example, let us consider a system of bodies to be 
oscillating about a position of equilibrium. We know by Lagrange's 
general solution, that the equation /(m) = is of an even order. 
Its roots are of the form 

Whatever the numerical values of these may be we have 

m^ + m^ + m^ = m^, m^+m^ + m^^ m^, &c. 

so that the small terms of the third, fifth, &c. orders might affect 
the stability of the oscillation. But we shall now show that they 
only affect the periods of oscillation, and not the stability of the 
system. 

Since both sides of Lagrange's equations must be of — 2 di- 
mensions in time and the impressed forces are also of — 2 dimen- 
sions, it is clear that these terms must consist of powers of x, y, &c, 

-j^^ j^> &c. and products of an even number of factors of 

dsc dij 

-ji* -j^y &c. We know also, by Lagrange's solution, that the 

co-ordinates take the form 

x= M^ cos {oit + XJ + &c. 
y = M^ cos {xt-\- \j) + &c. 
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d 

The minors a, a\ &c. are also all even functions of -r: • hence 

at 

the equation found after elimination is of the form 

/ ^^j « = P cos (at + \) + ^cos (fit + fi) + ... 

Replacing aV^, — aV— 1, &c. by m^, m,, &c. we find by 
Art. 6, 

P P 



Since /'(m) is of odd dimensions, and m^=:^m^ we clearly have 
Bm^=—8m^ and therefore the test is satisfied. 

9. [As another example let us apply the rule of Art. 6 to 
some very simple case which will involve no algebraical substi- 
tutions of any length. 

The motion of a simple pendulum under the action of gravity 
may be made to depend on the equation 

• g + aV = ;8a^ (1), 

where a and /8 are two constants and x is the inclination of the 
pendulum to the vertical which is supposed to be small. The first 
approximation to the motion is 

X == Mjfi^ + M^i^^ (2), 

where m^ and m^ are the roots of the equation m* + a' = 0. Our 
object is to ascertain by help of the rule given in Art. 6 whether 
the small force represented by Ba^ renders this first approximation 
unstable or merely slightly alters the numerical values of m^ and 



fn^* 



The two roots are connected by the relation 

^1 + ^8 = (3). 

Substituting the value of x on the right-hand side^ we have 

ide + *") ^ == ^^ ^i^^ [M^e'^^^ + -a/,e^*) + &c. 
Hence by the rule in Art. 6 

**^» 2m, Jf, ' ^* 2m,iJ/, ' 

These clearly satisfy the relation 

8m, + 8m, = 0, 
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and therefore the first approximation taken above, so far as the 
disturbing force Ba? is concerned, is stable. 

If the small force had been ^{-tl) instead of ^a?, it is easy 
to see in the same way that 

so that the relation Sm^ + S«ij=0 would not have been satisfied. 
The first approximation taken above is therefore not suflBciently 
accurate to serve as a first approximation throughout the motion. 

In this example we have considered the effect of a small force 
of the third order in disturbing the stability of the motion given 
by equation (!)• The same equation will obviously occur in many 
other cases of motion. For example, let a particle describe a cir- 
cular orbit about a centre of force situated in the centre. If 
slightly disturbed the equation giving the disturbance x in the 
radius vector takes the form 

where a, fi and 7 are constants. Similar remarks wiE therefore 
apply to this case also.] 

10. When the coefficients of the equation of motion are not 
strictly constant, but yet do not vary much, then we may transpose 
the small variable parts of these terms to the right-hand side of 
the equation, and treat their products by the differential coeffi- 
cients of the co-ordinates as small quantities of the second order. 
Suppose the variable part of one of these coefficients to be p sin nt, 
where p is small, and let/(m)=0 be the equation givmg the 
periods of oscillation of the system when the coefficients are taken 
constant. Then it is clear that unless n is nearly equal to the 
sum or difference of two values of m, this term cannot rise into 
importance. On proceeding to higher orders we see that these 
terms cannot produce important effects unless some commen- 
surable relation between m and the roots of the equation /(m)' = 
should be very nearly satisfied. 

11. It should be remarked that when the coefficients are not 
constant it is not a sufficient test of stability that they should 
always satisfy the conditions of stability obtained by giving them 
their instantaneotis constant values. Thus if the equation of 
motion were 
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the coefficient of x is always positive, yet as the equation is satis- 
fied hj x = a/^t,x may become as great as we please. 

Even if the coefficients are nearly constant, we must yet ex- 
amine, by the rules just given, if their small changes are so timed 
as gradually to increase the oscillation until the divergence from 
the given motion is no longer small. 

[Suppose a system to be oscillating so that its motion is de- 
termined by the equation 

^+20^ = 0, 

where j is a known function of t, which during the time under 
consideration always lies between ^ and ff* the latter being the 
greater. Let the system be started with an initial co-ordinate x^ 
and an initial velocity x^ in a direction tending to increase x. 
It may be shown that the system will begin to return, i.e. x will 

/ ^ 
begin to decrease before x becomes as great sis x/ x^ + -^. If 

± m, Tw' be two successive maximum values of x, we may also 

show that m! cannot be so great as ^ rn, and that the time from 

one maximum to the next lies between -5 and ^]. 



R. A. 



CHAPTER VIII. 



The HomUUmian Characteristic or Principal functions whenfonind deter- 
rminJe ai once the motion of the system from, one given position to 
another, and whether the motion is stable or unstable. Arts, 1 — 3. 

Examples with a made of effecting the integration S= \Ldt in small oscil- 
lations. Arts. 4 — 7, 

The Characteristic function supplies the condition tliat the motion is stable 
as to space only, while the Principal function gives the conditions that 
it is stable both as to space and time. Art, 8. 

In what sense the motion is unstable if either of the two ffamUtonicm 
functions is a minimum. Arts, 9 — li, 

1. If we had any convenient methods of finding the Hamil- 
tonian Characteristic or Principal function, we might determine 
without diflSculty the conditions of stability of a dynamical system 
at the same time that we deduce the integrals of the equations of 
motion. But it is very diflScult to discover either of these func- 
tions by an d priori method. We have indeed differential equa- 
tions which they must satisfy, and Jacobi has taught us what kind 
of solution will serve our purpose. But the difficulty of finding 
these solutions is as great as that of solving the equations of 
motion. For these reasons it does not seem necessary to dwell on 
the uses of these functions, 

2. Suppose the Principal function 8 of a, djmamical system to 
have been found in terms 5f the initial co-ordinates 0^ ^^ and the 
co-ordinates 0, (b and the time t. Let the semi-vis viva be given 

by 

T=P^-^ + Qd'<t>' + R^, 
where P, Q, R are known functions of 0, j). Let P^, Q^, R^ be 
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the values of these when 5^, <^^ are written for 6, ^. The final 
integrals of the equations of motion are then given by 






Let the ^tem receive any disturbance at the time < = 0, so 
that while starting from the same initial position, its initial ve- 
locities are slightly altered. Let x^, y^ be these initial chahges 
of 6^ and ^/ and let d + x, (f>+y be the co-ordinates of tho 
system at the time t. Then ure have 



Here x^, y^ are small arbitrary quantities, hence x and y will 
be ,..U if none »t tte r^io. of ^, ^. ^, or ^ 
to the determinant 

d^S tPS d^S d^8 



dd^dJd d^^d^ de^d^ d(f>^d0 ' 
be large. 

If the initial position as well as the initial motion be altered 
we may find, by a precisely similar process, the conditions that 
X and y should be small. If the system have more than two inde- 
pendent motions, we have more than two co-ordinates, but the 
conditions of stability are found in the same way. 

3. If 0? and y be small throughout the whole motion from the 
one given position to the other, not only does the system not 
deviate far from its undisturbed course, but the system at any 
instant is also very nearly coincident with its undisturbed place 
at the same time. It is important to notice this, for the word 
** stability " is sometimes used m a different sense. 

• 

This condition of stability may be put under a form in which 

no reference is made to 8. Let u and v be the components of 

momentum of the system corresponding to the co-ordinates d, d> 

dT dT 
respectively, i,e. let ** = 3^>> ^~^* ^^^ ^^^ these be expressed 
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as functions of Oq, <l>^f 6y 4> ^^^ ^- Then tlie preceding equation 
may be written 

du dv 

du dv ^ 

where a and fi are small arbitrary quantities. The condition of 
stability is that the values of x and y thus found should be small. 

4. As an example let us consider the case of a projectile. 
If 6 be the horizontal, and ^ the vertical co-ordinate of the par- 
ticle, we have 

The equations to find x and y are evidently 

1 



1 

Hence the system continually deviates more and more from its 
undisturbed place. 

5. In order to calculate the form of 8 when a system- is oscil- 
lating about a state of motion, it is convenient to choose as co- 
ordinates some small quantities x, y which vanish in the given 
state of motion. Let the Lagrangian function be written in the 
form 

where L^ is a homogeneous function of a?, y, x'y j/. Then by a 
theorem of Euler's 

where the 2*s imply summation for all co-ordinates. 

As in Art. 9 of Chap. IV. we have 

d dL^ dL^ 
dt dx dx * 

and the oscillations arc given by 

d dL,^ __ dL,^ 

dt did dx ' - 
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Hence we find 

Integrating we have 

Thus the integration has been effected, but in order to express 
S as a function of (c^f y^, x, y and t, it will be necessary to find 
a?/, y^', X, y in terms of these quantities. 

6. As an example, let the position of the system depend on 
one co-ordinate x and let 

where the coefficients are all constants. We then find by the 
process just indicated that il^ = and 

2 e^« — ^-«K * 

A 
where w' = -^ . Applying the criterion of stability we find that 

-1 — 7- will finally become small if m is real. The motion is there- 
axflx 

fore unstable or stable according as -4,i, B^ have the same or 

opposite signs. 

7. If the position of the systerfi depend on two co-ordinates 
a?, y, let 

L^L^-^-A^x-^Ajf-^-B^ai ^B^y' 

+ C^xx+C^t/ + G^,yx' + a^y'. 

We then find 
8=Ljt + B,{x'-x,) + B,{y-y,) 

where 
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If we now express a?', y\ a?/, y^ in terms of a?, y, x^ y^ and t, we 
find for a a fraction whose numerator is a homogeneous quadratic 
function of a?, y, x^, y^, the coeflBcients being linear functions of 
exponentials of ^, and whose denominator is another linear function 
of the same exponentials. These exponentials become sines and 
cosines when the motion is stable. Thus when the given motion 
is steady the simplest inspection of the form of <r will determine 
whether the motion is stable or not. 

Referring the motion to principal co-ordinates for the sake of 
brevity, and writing 2(y= (7^,— G^^ we find that <r must satisfy, 
the differential equation 

This equation is obviously satisfied by such a function as that 
just described. The solution of this equation may be reduced to 
linear equations and thus <r may be found. But it is unnecessary 
to dwell on this, for this would be equivalent to returning to the 
Hamiltonian equations. 

8. If we wish to determine the condition that the genieral 
course of a dynamical system is stable without requiring it should 
be near its undisturbed place at any the same time, it is more 
convenient to use the Characteristic function. Suppose that the 
Characteristic function has in Jacobi's manner been expressed as 
a function of the co-ordinates 0, ^, the constant A of vis viva and 
two arbitrary constants a^, a^. Then 

V=f{0,<l>,h,a,) + a,. 

The relation between and ^, which may be called the equa- 
tion to the path of the system, is given by 

dV , 

where \ is another constant. Let the system be disturbed from 
the same initial position so that the whole energy is unaltered. 
The change in corresponding to any given value of 5 is found 
from 

da^ 
d^V 



Let A be the initial value of -5— § , then 



^0 = — ^TpT'-^^i' 
da^d<f> 
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The condition that the path should be stable is that the coeffi- 
cient of Sttj should not be large. 

We might also use the function called Q by Sir W. R. Ha- 
milton, but it seems unnecessary to dwell more on this subject. 

9. The instability of a system may be deduced from the 
Hamiltonian Principal or Characteristic functions, expressed as a 
minimum. Suppose a dynamical system to move from one posi- 
tion A to another 5 in a time t, then the motion may be found by 

making the first variation o{ 8= j Ldt equal to zero, the time of 

transit being constant. The constants of integration are deter- 
mined by the conditions that the co-ordinates have given values 
when ^ = and t = t To determine whether ;S is a maximum or 
minimum or neither we must examine the second variation and 
here we have the assistance of Jacobi*s rule. The determination 
of the constants will depend on the solution of equations and may 
lead to several different kinds of motion from A to B, One of 
these will be the actual motion. Let us move B along this until 
one of the other motions coincides or as we may say approaches 
indefinitely near to this actual motion. We have then reached 
a boundary beyond which the integration must not extend if /S» is 
to be a maximum or minimum. See Todhunter's History of the 

cPL 

Calculus of Variations, page 251, Further ^^, if 5 be a co- 
ordinate, is positive throughout the limits of integration, so that 
S will be a minimum and not a maximum. 

10. When there are several co-ordinates 0, <^, &c, which are 
to be found as functions of the time, we may easily show that 
Jacobi's condition is a necessary one, and this is all that we require 
for the next proposition. If the system can move in two ways 
from A to B, then B8=^0 along each, and therefore when these 
two are adjacent we have both SS = and 8{8-hBS) = 0. This 
shows that the second variation can be made to vanish by taking 
one variation through the other. This second variation will then 
be the same as the quadratic term of the series obtained by 
changing the co-ordinates 0, <f> into + 80, <f>-{- S<l>, because we can 
take 8^0 = and S^^ = 0. Hence as the sum of the terms of the 
third order does not, in general, vanish for this displacement, it is 
clear that S cannot be either a maximum or a minimum. 

Let the actual motion be from A to B, and let a neighbouring 
motion starting from A lead the system to a position C reached in 
the same time along the actual motion before reaching B. Then 
we can show that a variation of the actual motion from A to B 
can be found which makes 8^8 of any sign. Let F be any positiou 



104 THE PRINCIPAL FUNCTION. [CHAP* 

on the neighbouring motion before reaching (7, and Q, one on the 
actual motion after passing G. Then considering P and Q as fixed 
and also the time of l^ansit^ the motion along PGQ cannot 
make 85= 0; for this condition is known to lead to the ordinary 
dynamical equations, and it is clear that (impulsive forces being 
set aside) no actual motion can be discontinuous. But there is 
discontinuity at (7, for otherwise when the system is started from 
B towards A, two courses would be open to the system on arriving 
at G. Hence the first variations of Sior an imaginary motion along 
PGQ are not zero, and therefore may have any sign. But since 
the discontinuity at G is of the first order of small quantities, this 
first variation is of the second order. Now the value of 8 for the 
actual motion is equal to that along the neighbouring motion to 
C and then along the actual motion to B. Hence, P and Q being 
still fixed, variations of the actual motion from ^ to J5 can be 
found which make B^8 of any sign. 

11. Let us apply this theory to determine the stability of a 
given state of motion. First let us suppose the given motion to 
be steady and to depend on only two variables, K we use the 
function S there will be one co-ordinate and the time, if V two 
co-ordinates. Let the system be disturbed at any moment by an 
alteration of the velocities of its several parts, so that the initial 
position of the disturbed motion is an undisturbed position. If 
the motion be stable the system will oscillate about the un- 
disturbed motion, the oscillation repeating itself at a constant 
interval. It follows therefore by Jacobi's rule that ;S or F cannot 
be a minimum for a period longer than the time of a half-oscilla- 
tion. If therefore /Sf or F be a minimum for all variations, starting 
from A and ending at B, where S is a position on the steady 
motion reached by the system at an interval as long as we please, 
then the motion is unstable. 

If we give a meaning to the word " stable" somewhat different* 
from its usual signification, we may extend this proposition to 
determine a test of the stability of any motion, whether steady or 
not. All we have assumed is, that, if the motion be not altogether 
unstable, there are some disturbances which will cause the system 
periodically to assume the same positions as it would have done if 
it had been undisturbed, but the interval of these periods may be 
any whatever provided the first be finite. If we use the Character- 
istic function, these disturbances must be such as not to alter the 
constant of vis viva, and if the Principal function, they must be 
such as to bring the system to an undisturbed position in the same 
time, 

* This meamng doea not always agree with the results of Art. 14, Chap, iv, 
[See also Arts. 17 and 18, Chap, iv.] ' 
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12. Next let us suppose that as the system proceeds from the 
initial position A along the actual motion, 8 ceases to be a 
miuimum at some position B. The conditions for a minimum 
are of two kinds. Suppose the system to depend on two co- 
ordinates dy <l>, and let L be the Lagrangian function, then (1) we 

must have -j^^ and -jffi^iii^'iffjTi ^oih positive, and (2) it 

must be possible to choose three arbitrary constants which enter 
into a very complicated expression, so that this expression may 
never become infinite between the limits of integration. The first 
condition is clearly always satisfied since the vis viva of any 
system is necessarily positive for all values of ff and ^'. The 
second condition will fail if there are two neighbouring motions by 
which the system can proceed from A to any position between 
A and B. If this be the mode of failure, it is clear from the 
reasoning of Art. 2 that the conditions of stability are satisfied for 
one kind of disturbance, and that therefore some at least of the 
harmonic motions are stable or oscillatory, though the motion may 
be unstable for a different kind of disturbance. 

13. [These conditions become much simpler when the position 
of the system is determined by one co-ordinate, or when the 
Lagrangian function can be reduced to depend on one co-ordinate. 
Let this co-ordinate be so chosen that it vanishes along the given 
motion, and let us also suppose that both it and its differential 
coefficient with regard to t^ are small for all neighbouring con- 
strained motions. Let this co-ordinate be called 6 and let the 
Lagrangian function be 

L^L, + A,e ^B^ff+\A,,0' + \ BJ'' + CJff. 

Then since the Lagrangian equation of motion is satisfied by 
hypothesis when ^ = 0, we have A^ = B( where the accent, as 
usual, denotes differentiation with regard to t 

If the system be now conducted from the initial position A to 
any other position B, both on the given motion, by any neighbour- 
ing mode of motion, we have 

Sf = JL,dt + J(i BJ' + CJff + \ AJP^ dt 
If ^ = w be any solution of the Lagrangian equation 

R. A. 8 
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we may write the function 8 in the form* 

The second term is essentially positive, since B^^ must be 
positive. Hence iSf is a minimum along the given motion unless 

we can so choose the arbitrary displacement as to make -^ - = 0. 

This gives d^cu where c is some constant. But 6 must vanish 
at the two limits A and B, hence this choice of 6 is excluded 
unless there is some neighbouring mode of motion by which the 
system could move freely from the given initial position A to the 
position B, The result is that S cannot cease to be a minimum 
before the first instant at which some neighbouring motion will 
bring the system (starting from A) into coincidence with some 
contemporaneous position on the given motion. If the given 
motion be steady, it follows that S cannot cease to be a minimum 
before a time which is half that of a complete oscillation. 

We thus have a test of stability. If the system depend on one 
co-ordinate and if ^ he a Tninimwm when the limits of integration 
are from the initial position A to all positions on the actual motion^ 
that motion is unstable. But if S cease to be a minimwm at some 
point C, then the actttal motion is stable from A to C] 

* [Following Lagrange's rule we may write the second term of iS in the form 

The quantity ontside the integral sign is to he taken hetween the given limits 
and is zero, since d yanishes at each limit. Let us now pat 

It is clear that this valne of X cannot he infinite hetween the limits of integra- 
tion nnless u vanishes. For hy hypothesis u an4 u' are hoth finite and the co- 
efficients B^i and C^ in the Lagrangian function are also finite. It then easily 
follows from the equation _ 

^^ ( jBj jU + Oiiu) = CiiU + A iiti, 

that {On + 2X)>= Bu {A^ + 2X'). 

Hence the second term of S hecomes 



lJBn{e--e^yst. 



which is the result in the text. This might also have heen deduced from Jacobi's 
general transformation with one independent variable gtven in Prof. Jellett's 
CalculiL8 of Variations or Prot Price's Differential Calculus, 

If u=0 between the limits of integration this transformation fails, but, as is 
evident from the argument in the text, we choose ^=w to represent a neighbouring 
free motion such that u=sO just before the system reaches A ; also B is so placed 
that the next instant at which u=0 is after the system has passed B.] 
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14. [We shall conclude the chapter with the application of 
this criterion of stability to some simple case, 

A particle describes a circular orbit about a centre of force 
situated in the centre. It is required to deduce the conditions of 
stability as to space from the Characteristic fum^tion V. 

Let a be the radius of the circle, n the angular velocity of the 
particle about the centre 0. Let ^(a)be the law of force. Let 
A and B be two points taken on the circular orbit, and let the 
particle be conducted from -4 to £ by some neighbouring path 
with the same energy as in the circular orbit. Let r=a + phe 
the radius vector of this path, corresponding to any angle 0, 

If V be the velocity at any point of this path, we easily find 

If s be the arc of the path, we have 

ds . ^ 1 fdpV 

If the angle A OB = 13, we therefore have 
r^jvds = a^n^^lf^{[^^^^^ 

where/ = ^+3. 

If the neighbouring path be a free path described with the 
same energy, its equation is 

p = L sin p0, 

where is measured from the radius vector OA, and L is an 
arbitrary constant. This free path will cut the circle again in 
some point C. J£ the angle -4 0(7= 7, we have py^ir. 

If B coincide with C, we find by substituting this value of p in 
the expression for V, that the second term of Fis zero. If ^ be 
beyond C but such that the angle GOB. is less than 7, draw two 
free paths one from A as before and the other backwards from B 
to meet the former in some point P. Then the angle A OP = ^13. 
If the particle be conducted from -4 to P aloDg one path and 
from P to P along the other, we find that the excess of the 
action over the action in the circular arc AB is equal to 

-—jL'sin^^. 

Since p^ is greater than ir and less than 2'rr, this excess is 
negative. The action along the circular arc is therefore not a 
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minimum if J& be beyond the first intersection of a neighbouring 
free path. 

Lastly we may show that the action is a minimum if J3 lie 
between A and G. To prove this we write the integral in the 
expression for V in the form 

The first of these two terms is zero since p vanishes at each 
limit Following Lagrange's rule we make 

The integral then becomes 

This is always positive and the action along the circular arc is 
a minimum. The argument however requires that X should not 
become infinite between the limits of integration. It is easy to 
see that 

where £ is a constant to be chosen at our pleasure. But if fi 

IT 

exceed - it is impossible to choose E bo as to keep X finite be- 
tween the limits of integration. Hence the action is a minimum 
only if the angle A OB subtended by the limiting positions at the 

centre is less than - . The circular orbit is therefore stable if »* 

be positive. 

If, however, p* be negative, the expression for X changes its 
character. Writing — j" for ^*, we find 

q — X ' 

thus the value of X can be chosen so as not to be infinite for all 
positive values of 0. In this case the function F is a minimum 
for all arcs AB however distant B may be from A. The circular 
motion is therefore unstable.] 
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